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Investigation  of  the  interference  heating  phenomenon  was  conducted  at  Mach  number  8 
over  a  Reynolds  number  range  from  .  4  x  10®  to  3. 7  x  10b  per  ft  In  Hypersonic  Wind  Tunnel 
B  of  the  von  Karman  Gas  Dynamics  Facility,  Arnold  Engineering  Development  Center#  Shock 
generator  models  consisted  of  1.5,  5,  10  and  15  deg  wedges;  2.5,  7.5  and  12.5  deg  cones; 

.  5  and  2. 0  in.  diameter  hemispheres  and  an  orbiter-type  vehicle.  Shock  receiver  models 
included  sharp  and  blunt  flat  plates  and  a  booster  vehicle.  Both  pressure  distributions  and 
heat  transfer  distributions  (thermocouple  and  phase  change  paint)  were  obtained  in  the 
regions  of  shock  impingement.  The  test  program  was  organized  in  a  building-block  approach. 
A  limited  number  of  sharp  plate  runs  were  performed  with  a  boundary  layer  trip  to  obtain 
turbulent  data  for  comparison  with  other  experiments.  Our  heating  amplification  for  two- 
dimensional  shocks  agreed  with  previous  results  that  ' 
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Turbulent  correlations  were  also  obtained  for  three-dimensional  conical  and  spherical 
shocks.  The  majority  of  our  testing  used  wedge  shocks  impinging  on  the  sharp  plate  with 
an  initially  laminar  boundary  layer.  A  wide  range  of  Reynolds  numbers  were  covered  by 
testing  at  four  tunnel  pressures  and  three  x  locations  on  the  plate.  Laminar  interactions 
correlated  as 
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ABSTRACT 


Investigation  of  the  interference  heating  phenomenon  was  conducted  at  Mach  number  8 
over  a  Reynolds  number  range  from  .  4  x  10°  to  3.  7  x  106  per  ft  In  Hypersonic  Wind  Tunnel 
B  of  the  von  Karman  Gas  Dynamics  Facility,  Arnold  Engineering  Development  Center.  Shock 
generator  models  consisted  of  1.5,  5,  10  and  15  deg  wedges;  2.5,  7.5  and  12.5  deg  cones; 

.  5  and  2.0  in.  diameter  hemispheres  and  an  orbiter-type  vehicle.  Shock  receiver  models 
included  sharp  and  blunt  flat  plates  and  a  booster  vehicle.  Both  pressure  distributions  and 
heat  transfer  distributions  (thermocouple  and  phase  change  paint)  were  obtained  in  the 
regions  of  shock  impingement.  The  test  program  was  organized  in  a  building-block  approach. 
A  limited  number  of  sharp  plate  runs  were  performed  with  a  boundary  layer  trip  to  obtain 
turbulent  data  for  comparison  with  other  experiments.  Our  heating  amplification  for  two- 
dimensional  shocks  agreed  with  previous  results  that 

hpeak  _  ( ppeak 
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Turbulent  correlations  were  also  obtained  for  three-dimensional  conical  and  spherical 
shocks.  The  majority  of  our  testing  used  wedge  shocks  impinging  on  the  sharpplate  with  an 
initially  laminar  boundary  layer.  A  wide  range  of  Reynolds  numbers  were  covered  by  test¬ 
ing  at  four  tunnel  pressures  and  three  x  locations  on  the  plate.  Laminar  interactions  corre¬ 
lated  as  / 
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Interactions  that  resulted  in  a  turbulent  boundary  layer  after  impingement  correlated  as 
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The  effects  of  a  very  blunt  leading  edge  and  three-dimensional  shock  generators  on  laminar 
interactions  were  studied  and  some  trends  noted,  but  the  data  are  not  extensive  enough  for 
definite  conclusions.  Finally,  the  major  correlations  and  trends  resulting  from  the  simple 
geometry  testing  were  applied  to  the  orbiter /booster  configuration  and  compared  to  the 
actual  test  data. 
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Section  I 


INTRODUCTION 

The  increase  in  heating  that  occurs  when  a  shock  wave  impinges  on  a  vehicle  surface 
can  be  a  major  design  problem  at  high  speeds.  Consequently,  the  subject  has  received 
increasing  attention  as  aircraft  and  missile  speeds  increase.  Ryan,  in  Ref.  1,  summarizes 
the  analytical  and  experimental  literature  from  1947  to  1968.  Shock  interaction  with  a 
turbulent  boundary  layer  has  been  treated  mostly  by  experiment.  Sayano,  in  1962  fRef  2). 
presented  an  empirical  correlation  of  the  heating  amplification  as  a  function  of  the  pressure 
amplification  which,  in  turn,  can  be  determined  by  a  straightforward  inviscid  calculation. 
Several  analytical  approaches  to  the  laminar  case  have  been  proposed.  Murphy,  in  Ref  3, 
reviewed  three  of  the  methods,  and  concluded  that  they  agreed  with  the  data  at  low  Mach 
numbers  and  shock  strengths  but  deviated  for  strong  shock  waves  at  hypersonic  speeds. 

Rose,  in  Ref  4,  proposed  a  two-layer  approach  applicable  to  laminar  or  turbulent  boundary 
layers.  While  ail  the  analytical  methods  provide  insight  into  the  problem,  they  are  not  yet 
refined  to  the  point  of  being  applicable  for  design  calculations  at  hypersonic  speeds. 
Markarian,  therefore  in  Ref  5,  presented  empirical  correlations  of  heating  amplification  data 
for  various  shock  strengths  and  laminar,  transitional  and  turbulent  boundary  layers.  The 
purpose  of  our  study  is  to  extend  the  empirical  approach  of  Sayano  and  Markarian  to  provide 
the  design  engineer  with  a  simple  method  of  making  preliminary  heating  estimates. 

Turbulent  amplifications  are  relatively  well  behaved  and  several  investigators  have  proposed 
the  correlation, 

kpeak 

hUnd 

where  N  is  .8  to .  85.  We  therefore  took  only  limited  turbulent  data,  primarily  to  show  consistency 
with  other  data,  although  some  new  data  with  three-dimensional  shocks  were  also  taken. 

The  main  emphasis  of  our  program  is  to  investigate  the  laminar/transitionai  regime  at  a 
hypersonic  Mach  number.  Our  goal  is  a  correlation  of  flat  plate  heating  amplification  as 
a  function  of  Reynolds  number  and  wedge  shock  strength.  Since  the  design  engineer  faces 
practical  problems  where  the  incident  shock  cannot  be  considered  two-dimensional,  we 
have  also  taken  data  using  cone  and  sphere  shock  generators.  The  effect  of  plate  pressure 
gradient  is  investigated  using  a  blunt  leading  edge.  Finally  we  will  treat  the  practical 
problem  of  two  airplane  shapes  (orbiter  and  booster)  in  close  proximity.  Predictions 
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Section  II 


EXPERIMENTAL  PROGRAM 

This  section  will  summarize  the  model  design  and  the  test  program.  Detailed 
discussion  and  drawings  of  the  models  are  contained  in  the  test  report  (Ref  6). 


1.  DESCRIPTION  OF  MODELS 

The  models  consist  of  a  set  of  stainless  steel  shock  generators  which  include  one 
6  x  9  in.  fiat  plate  wedge,  two  hemisphere-cylinders  having  respective  diameters  of  2. 0  and 
.  5  in. ,  one  2. 5  deg  semi-vertex  angle  cone-cylinder  with  a  .  5  in.  base  diameter,  two  7. 5  deg 
semi  vertex  angle  cone-cylinders  with  respective  base  diameters  of  1. 5  in.  and  2. 5  in. ,  one 
12. 5  deg  semi-vertex  angle  cone-cylinder  with  a  2. 5  in.  base  diameter,  and  an  aerodynamic 
configuration  representative  of  a  typical  shuttle  orbiter  vehicle  having  a  spherical  nose 
radius  of  .25  in.  These  generators  are  shown  schematically  in  Fig.  3  thru  7. 


Each  generator  may  be  positioned  above  any  one  of  three  instrumented  shock  receiver 
configurations:  namely,  a  flat  plate  having  interchangeable  blunt  and  sharp  leading  edges; 
a  1.0  in.  diameter  hemisphere-cylinder  (not  used  in  the  current  program);  and  an  aero¬ 
dynamic  configuration  representative  of  a  typical  shuttle  booster  vehicle  (See  Fig.  8  and  9). 

Schematics  of  typical  generator-receiver  arrangements  are  shown  in  Figures  10  and  11. 
All  support  members  are  fabricated  from  stainless  steel.  Wedge  generator  angles  of  attack 
were  selected  at  1. 5,  5,  10  and  15  deg  by  rotating  the  support  blade  sector  with  respect  to 
the  support  blade.  The  sector  rotates  about  a  fixed  point  relative  to  the  support  blade,  which 
is  coincident  with  the  reference  point  of  each  shock  generator  when  that  generator  is 
mounted  on  the  sector. 

The  model  support  system  allows  independent  horizontal  and  vertical  positioning  of 
the  generators  above  the  receiver.  The  support  blade  is  free  to  slide  vertically  in  the 
forward  end  of  the  horizontal  support  arm.  A  total  of  7  inches  of  vertical  travel  in  0. 1-in. 
increments  is  possible.  The  shaft  of  the  support  arm  is  free  to  slide  horizontally  in  a  split 
clamp  located  at  the  upper  end  of  the  support  stand,  allowing  a  total  of  20  in.  of  horizontal 
travel  in  0. 1-in.  increments. 

The  flat  plate  shock  receiver  (Fig.  8)  is  essentially  a  frame  designed  to  accept  three 
interchangeable  24  x  24  x  .  5  in.  instrumented  inserts  for  the  separate  acquisition  of  pressure. 


3 


‘l 

% 


nliHiT'l . 


thin-skin  heat  transfer  and  phase  change  paint  heat  transfer  distribution  data.  The  receiver 
frame  will  also  accept  either  a  blunt  or  sharp  leading  edge. 

Typical  cross-sections  of  the  thin  skin  (thermocouple)  and  pressure  inserts  are  shown 
in  Fig.  12.  The  thin-skin  insert  consists  of  a  nominal  .  050-in.  thick  sheet  of  AISI  321 
stainless  steel  mounted  on  a  slab  of  Teflon.  Low  friction  slide  washers  and  lock  nuts  are 
used  to  fasten  the  skin  to  the  slab  as  shown  in  Fig.  12a.  This  permits  the  skin  to  slide  on 
the  Teflon,  which  is  free  to  slide  on  the  frame.  Such  a  design  was  utilized  for  expansion 
compensation.  The  stainless  steel  sheet  is  instrumented  with  121  laser  welded  thermo¬ 
couples  primarily  along  its  longitudinal  centerline.  A  .  75-in.  wide  channel  is  cut  through 
the  Teflon  slab  in  instrumented  areas  to  minimize  heat  conduction  from  the  skin  into  the 
Teflon  slab  in  this  region. 

The  pressure  insert  is  a  solid  stainless  steel  plate  instrumented  with  121  pressure 
taps  at  the  same  locations  as  the  thermocouples  of  the  thin-skin  insert.  The  insert  is 
bolted  to  the  receiver  frame  as  shown  in  Fig.  12b.  Slide  washers  are  used  to  provide 
compensation  for  thermal  expansion. 

The  phase  change  paint  insert  is  a  solid  slab  of  black  Teflon  which  is  secured  to  the 
frame  in  the  same  manner  as  the  thin  skin  insert.  Ten  thermocouples,  five  to  either  side 
of  the  longitudinal  centerline,  are  also  embedded  in  the  slab  for  initial  temperature  measure¬ 
ment. 

All  three  flat  plate  inserts  are  equipped  for  installation  of  a  boundary  layer  trip 
device  1. 5  in.  from  their  forward  edges.  This  locates  the  trip  4.0  in.  aft  of  the  model 
leading  edge  when  the  inserts  are  installed  on  the  receiver  frame  (See  Fig.  8).  The  trip 
is  constructed  of  .125-in.  diameter  ball  bearings  set .  371 -in.  apart  from  center  to  center, 
which  are  recessed  and  soldered  into  beveled  steel  strips.  Due  to  the  recess,  the  effective 
trip  height  with  respect  to  the  strip  upper  surface  is  .  081  in.  The  trips  are  installed  for 
selected  sharp  flat  plate  runa  to  investigate  shock  interaction  with  an  initially  turbulent 
boundary  layer. 

Two  booster  vehicle  shock  receivers  of  identical  configuration  are  available.  The 
first  is  fabricated  from  black  Teflon.  This  model  is  used  solely  for  the  acquisition  of 
phase  change  paint  heat  transfer  distribution  data.  A  reference  grid  pattern  of  embedded 
white  Teflon  rod  is  incorporated  to  facilitate  data  reduction  (See  Fig.  13).  The  second 
booster  model  is  fabricated  from  stainless  steel  and  is  designed  for  the  separate  acquisition 
of  pressure  and  thin-skin  heat  transfer  distribution  data.  The  model  consists  of  three 
sections,  i.e. ,  nose,  body  and  wing.  (See  Fig.  9).  Three  interchangeable  nose  pieces  are 
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available:  the  first  is  pressure-instrumented  and  contains  56  taps;  the  second  is  instru¬ 
mented  with  118  thermocouples,  and  the  third  is  a  spare  containing  no  instrumentation 
at  all.  A  single  wing  section  has  provision  to  accept  interchangeable  inserts  as  indicated 
by  the  cross-hatched  areas  of  Fig.  9.  Three  are  provided:  a  pressure  insert  with  44  taps; 
a  thin-skin  insert  with  82  thermocouples;  and  an  uninstrumented  spare.  The  instrumented 
wing  inserts  are  similar  in  design  to  those  of  the  flat  plate  receiver.  The  booster  wing 
and  nose  sections  attach  to  a  common  body  section  which  accepts  the  support  sting. 


Two  crbiter  vehicle  shock  generators  of  identical  configuration  are  also  available. 

The  first  is  fabricated  from  black  Teflon.  This  model  is  tested  with  the  Teflon  booster 
receiver  for  the  determination  of  mutual  interference  heat  transfer  distributions  using  the 
phase  change  paint  technique.  The  model  incorporates  an  embedded  white  Teflon  reference 
grid  pattern  as  shown  in  Fig.  14.  The  second  orbiter  model  is  fabricated  from  stainless 
steel  and  is  designed  for  the  separate  acquisition  of  pressure  and  thin-skin  heat  transfer 
distribution  data.  This  model  may  also  be  used  as  an  uninstrumented  shock  generator. 

As  in  the  case  of  the  booster  wing  section,  the  steel  orbiter  can  accept  three  interchangeable 
inserts:  pressure  instrumented,  thermocouple  instrumented  and  uninstrumented  spare. 

The  pressure  insert  contains  100  taps  and  the  thin  skin  insert  200  thermocouples.  Insert 
location  is  indicated  by  the  cross-hatched  areas  of  Fig.  7. 

a.  Instrumentation 

Thirty  gage  chromel-alumel  premium  grade  wire  was  used  for  all  thermocouple 
instrumentation.  All  thermocouples  used  to  acquire  thin-skin  data  were  laser  welded  to 
nominal  .050  in.  thick  AIS1  321  stainless  steel. 

The  thin-skin  flat  plat*1  insert  contains  121  thermocouples  located  and  identified  as 
shown  in  Fig.  15.  One  hundred  sixteen  of  these  are  iD  tailed  .2  in.  apart  along  the  longi¬ 
tudinal  centerline  of  the  insert.  The  remaining  five  are  spaced  .  8  in  apart  along  a  line 
normal  to  the  centerline  at  a  station  21.5  in.  aft  of  the  insert's  front  edge  (24  in.  from  the 
leading  edge). 

Ten  thermocouples  are  embedded  in  the  Teflon  phase  change  paint  flat  plate  insert. 
These  are  used  for  initial  temperature  measurement  of  the  material  prior  to  model  injection. 
Five  thermocouples  are  positioned  5  in.  apart  along  longitudinal  lines  located  2.5  in.  to 
either  side  of  the  insert  centerline. 

The  thin-skin  orbiter  insert  contains  200  thermocouples  located  and  identified  as 
shown  in  Fig.  16.  They  are  positioned  along  six  rays  which  emanate  from  a  point  on  the 


5 


V 


r* **t&  -v?;r-.  ■  rw-^^v'-^r — -“-•sr-?—*:-— . jv  — — tt— 5^- r"^?^ ~ - - 


longitudinal  centerline  located  2.253  in.  aft  of  the  nose  section.  Thermocouples  are 
spaced  .2  in.  apart  along  their  respective  rays. 

The  locations  of  booster  thermocouples  and  their  numerical  identification  are  shown 
in  Fig.  17.  The  instrumented  nose  piece  contains  118  thermocouples  which  are  distributed 
.2  in.  apart  along  three  longitudinal  rays  and  10  deg  apart  along  the  circumferences  of  two 
cross-sections.  The  longitudinal  rays  emanate  from  the  nose  origin  and  are  defined  by 
the  intersection  of  the  booster  lower  surface  contour  and  the  0=0,  60,  and  90  deg  radial 
planes.  Instrumented  cross-sections  are  located  1. 983  in.  and  7.957  in.  aft  of  the  nose 
origin. 

The  booster  wing  thin -skin  insert  is  instrumented  with  82  thermocouples  as  shown  in 
Fig.  17.  Thirty-four  of  these  lie  along  each  extension  of  the  0  =60  deg  and  6  =90  deg  nose 
piece  longitudinal  rays.  The  remaining  14  lie  along  a  line  normal  to  these  extensions  at  a 
station  15. 6  in.  aft  of  the  booster  nose  origin. 

Stainless  steel  tubing  of  .  063  in.  outside  diameter  and  .  012  in.  wall  thickness  was 
used  for  all  pressure  instrumentation.  This  tubing  is  silver  soldered  to  larger  tubing  of 
,  093  in.  outside  diameter  and  .  012  in.  wall  thickness  approximately  1  ft.  aft  of  the  base 
of  the  assembled  models.  The  extension  tubing  is  18  ft.  in  length. 

The  flat  plate  pressure  insert  contains  121  pressure  taps  whose  locations  are  identi¬ 
cal  to  those  of  the  thermocouples  installed  in  the  thin-skin  insert.  Additional  pressure 
taps  are  installed  in  the  interchangeable  blunt  and  sharp  flat  plate  receiver  nose  pieces. 

Each  piece  contains  two  pressure  taps  located  1  in.  and  2  in.  aft  of  the  leading  edge  stations. 
This  instrumentation  is  in  line  with  the  insert  instrumentation  when  the  nose  pieces  are 
installed  on  the  receiver  frame  (See  Fig.  15). 

One  hundred  pressure  taps  are  installed  on  the  orbiter  pressure  insert.  These  are 
distributed  along  the  longitudinal  centerline  (Ray  3)  and  Rays  1  and  2  (See  Fig.  16),  and 
duplicate  the  thermocouple  locations  on  those  ;*ays. 

The  pressure-  instrumented  booster  nose  contains  56  taps  which  are  distributed  along 
longitudinal  rays  determined  by  the  9=90  and  120  deg  radial  pianes.  Forty-six  of  these 
are  installed  at  locations  identical  to  those  of  thermocouples  324-369  on  the  thin-skin 
booster  nose.  The  remaining  10  are  uniquely  positioned  .  4  in.  apart,  symmetrically 
opposite  to  the  locations  of  thermocouples  431,  433,  435,  437,  439,  441,  443,  445,  447 
and  449  (See  Fig.  17). 

The  booster  wing  pressure  insert  contains  44  taps  which  are  positioned  along  the 
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extensions  of  the 5=90°  and  120°  nose  piece  longitudinal  rays.  Thirty-four  of  these 
are  installed  at  locations  identical  to  those  of  thermocouples  370-403  on  the  thin-skin 
booster  wing  insert.  The  remaining  ten  are  uniquely  positioned  .  4  in.  apart  symmetrically 
opposite  to  the  locations  of  thermocouples  4")!,  453,  455,  457,  459,  461,  463,  455,  467  and 
469  as  shown  in  Fig.  17. 


Model  surface  pressures  were  measured  by  1  psid  and  15  psid  trarsducers. 


Phase  change  paint  data  pictures  were  obtained  with  70  mm  Varitron  cameras  at 
2  frames  per  second  framing  rate  up  to  15  seconds;  after  15  seconds  the  rate  was  reduced 
to  1  frame  per  second.  Fluorescent  bulb  light  banks  were  used  to  ilhiminate  the  test  section. 
The  film  used  was  Kodak  TR1-X  Pan  black  and  white.  The  time  of  each  picture  taken  was 
recorded  on  magnetic  tape  and,  when  correlated  with  wind  tunnel  conditions,  was  used  to 
calculate  the  heat-transfer  coefficient  corresponding  to  the  particular  elapsed  time  and  melt 
temperature  of  the  paint.  Fig.  18  shows  the  camera  setups  used  to  acquire  phase  change 
paint  data. 


2.  TEST  CONDITIONS 


Prior  to  each  run,  the  shock  generator  was  set  at  the  proper  angle  of  attack  and 

positioned  above  the  receiver  at  a  predetermined  location.  Generator  locations  were 

defined  by  the  horizontal  and  vertical  coordinates  of  the  generator  reference  point  denoted 

as  Xz  and  Z  ,  respectively.  Reference  coordinate  systems  applying  to  each  receiver 
m  m 

are  shown  in  Fig.  19  and  20,  The  models  were  injected  into  the  flow  following  the  estab- 
;  lishment  of  the  required  test  conditions.  Injection  time  from  fully  retracted  to  fully  injected 

positions  was  2.05  seconds.  The  approximate  locations  in  the  test  section  of  the  flat  plate 
i  and  booster  receivers  during  data  acquisition  are  shown  inFig.  19  and  20.  Receiver  align¬ 

ment  with  the  flow  was  monitored  visually  by  use  of  an  externally  located  scope.  Adjust¬ 
ments  were  accomplirhed  by  varying  the  angle  of  the  main  sting  support  system.  Alignment 
was  adjusted,  following  tunnel  pressure  changes,  to  compensate  for  attitude  variation  due  to 
t  changing  air  loads. 

All  thin-skinned  model  installations  were  followed  by  a  continuity  check  of  the  individual 
thermocouples.  A  thin-skin  run  was  initiated  by  model  injection  and  the  data  acquisition 
r  period  varied  from  5  to  10  seconds  depending  on  the  configuration.  Shadowgraphs  were 

■'  taken  through  either  or  both  the  viewing  ports  while  data  was  being  recorded.  Following 

?  model  retraction  from  the  test  section,  high  pressure  air  was  used  for  cooling. 


7 


All  pressure  model  installations  were  followed  by  leak  cheek  of  the  individual  pressure 
taps.  A  pressure  run  was  initiated  by  model  injection.  The  model  remained  in  the  test 
section  while  stabilization  of  the  surface  pressures  were  monitored  in  the  control  room. 
Stabilization  periods  were  lengthy  and  varied  according  to  the  configuration.  Once  the 
pressures  were  considered  stabilized,  data  acquisition  was  initiated.  Schlieren  photographs 
were  token  through  either  or  both  of  the  viewing  ports  while  data  was  being  recorded.  A 
phase  change  paint  run  was  initiated  by  model  injection.  Cameras  were  activated  just 
prior  to  the  time  the  model  reached  its  fully  injected  position  in  the  test  section.  Data  was 
acquired  over  about  a  60  second  interval  and  filming  was  terminated  as  model  retraction 
began.  High  pressure  air  was  used  for  cooling  following  each  run. 

In  all,  a  total  of  290  thermocouple,  lxl  pressure  and  41  phase  change  paint  runs  were 
made.  Table  I  shows  an  overview  of  the  run  schedule  and  the  four  nominal  free  stream 
test  conditions. 

3.  AEDC  DAT  '.  REDUCTION 

a.  Thermocouple  Heat  Transfer  Data 

The  reduction  of  thin  skin  thermocouple  data  normally  involves  only  the  calorimetric 
heat  balance  which,  in  coefficient  form,  is 

"d  Tw/df 

hT  =  *3TLBCp,STL  Tq-Tw  W 

v  «ar 

2  O 

hT=  Heat  transfer  coefficient  BTU/ft  -hr-  R 

p  STL=  Densi(y  of  stainless  steel  model  skin,  lb/ft3 
B  =  Model  skin  thickness  at  thermocouple  location,  ft 
Cp  STL=  Specific  heat  of  stainless  steel,  BTU/lb-  °R 
dTw/dt  =  slope  of  temperature-time  curve.  °R/  hr. 
t  -  time,  hr. 

Tq  =  Tunnel  Total  Temperature,  °R 
Tw  =  Skin  temperature,  °R 

Radiation  and  conduction  losses  are  neglected  in  this  heat  balance  and  data  reduction 
simply  requires  evaluation  of  dTw/dt  from  the  tsmperature-time  data  and  determination  of 
model  material  properties.  For  these  tests,  radiation  effects  were  negligible;  however, 
lateral  conduction  in  the  thin  skin  was  potentially  significant  in  the  impingement  regions  of 
the  models  where  large  gradients  in  heating  rates  were  expected. 
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C. ti! ^  -^" 


Separation  of  variables  and  integration  of  Equation  (1)  assuming  constant  f>  grpL» 


B*  Cp,  STL  md  To  -rields 


fi  B  C  i  ^“V  =  -^n  Pf 

STL  *  p,  STL  |_T 

where  i  refers  to  initial  conditions. 


T  -  T 

O  W,i| 

r^r — 

o  W  J 


Differentiation  of  Equation  (2)  with  respect  to  time  gives 


=  1 A 


T  -  T  . 
o  w,i 

T  -  T 
o  w 


STL  p,  STL  L  \  °  w  /  J  rT  _T  T 

Since  the  left  side  of  Equation  (3)  is  a  constant,  plotting  Jtn  —2 — versus  time 
will  give  a  straight  line  if  conduction  is  negligible.  Jo-Tw  _ 


IT  -T 

L  o  w 


The  data  were  evaluated  in  this  manner  and,  generally,  a  linear  portion  of  the 
curve  was  found  for  all  thermocouples.  For  high  heating  rates,  such  as  experienced 
in  the  impingement  regions,  the  linear  portion  was  quite  short.  A  linear  least  squares 
curve  fit  of  Jn  versus  time  was  applied  to  the  data  beginning  at  the  time 


curve  fit  of  Jn  W3-  ■  versus  time  was  applied  to  the  data  beginning  at  the  time 

^o"  w 

at  which  the  model  reached  the  tunnel  centerline  and  extending  for  a  time  span  which 
was  a  function  of  the  heating  rate,  shown  below: 

Heating  Rate,  °R/sec  Time  Span  of  Data  Number  of  Data 


Time  Span  of  Data 
_ Used.  Sec 


Number  of  Data 
Points  Used 


4<dTw/dt  0.5  11 

2<dTw/dt  <4  0.6  13 

dTw/dt  <2  1.0  21 

dTw/dt  <1  2.0  41 

Strictly,  the  value  o'  CP  STL  is  not  constant  as  assumed,  and  the  relation 

°p,  STL  “  Ao  +  A1  <  V  +  A2  < V2  +  A3  < Vs-  <4> 

was  used  with  appropriate  coefficients  and  the  value  of  Tw  at  the  midpoint  of  the  curve 
fit.  The  maximum  variation  of  C  OTT  over  any  curve  fit  was  less  than  one  percent;  thus 

p.  Oil.  g 

the  assumption  of  constancy  was  not  grossly  violated.  A  constant  494  lb^/ft  was  used 
for  Pg-pL  anc*  measured  values  of  skin  thickness,  B,  for  each  thermocouple  were  used. 
Use  of  the  data  reduction  equations  in  coefficient  form,  and  restriction  of  the  fit  times  to 
linear  portions  of  the  curve  should  preclude  the  necessity  for  correcting  for  lateral  con¬ 
duction  errors. 


\  J 


The  data  include  a  wide  range  of  heat-transfer  rates,  the  lowest  of  which  have  tem¬ 
perature-time  slopes  which  are  at  the  lower  limit  of  descrimination  of  the  data  acquisition 
system  (0.2  deg  R/sec).  The  uncertainty  of  the  data  is  a  function  of  the  temperature  slope, 
and  Fig.  21  illustrates  the  uncertainties  encountered.  First  the  uncertainty  (lcr )  of  the 
straight-line  data  fit  is  shown.  Also  shown  is  the  uncertainty  curve  for  the  total  system, 
including  test  condition  variations,  variations  in  model  pitch  angles  caused  by  wind  loads 
for  different  configurations,  and  plate  surface  irregularities  occurring  as  a  result  of 
shock  impingement. 

b.  Phase  Change  Paint  Keat  Transfer  Data 

Phase  change  paint  data  were  reduced  according  to  the  semi-infinite  slab  solution 
of  the  transient  one-dimensiona-  heat  conduction  equation  as  developed  in  Ref.  7.  The 
solution  is  given  as 
B2 

T  =  1  -  e  p  erfc  0 


Where  the  complementary  error  function  is  given  as 


erfc/? 


=1  P 


OO  2 


The  temperature  parameter  T  is  defined  as 

T  -T. 

—  _  pc  i 

1  "R/T  }-T  ‘ 
o'  i 

where  i  refers  to  initial  conditions  and 

T  =  phase  change  temperature  of  the  paint,  °R 
pc 

R  =  Recovery  factor 

The  parameter,  T ,  therefore,  is  a  constant  for  a  particular  run.  Corresponding  to 
each  value  of  T  there  exists  a  singular  value  of  0 

This  parameter  was  used  to  compute  the  heat  transfer  coefficient  corresponding  to  a 
phase  change  line  at  any  time  during  the  run  from  the  equation 
F/p  \  /r  \  FI  .  5 


,[(PTEf)  (Cp,  TEf)  (kTEF^j 


Where 


h  =  Heat  transfer  coefficient,  BTU/ft2-  hr-°R 

3 

Prrw  =  Density  of  Teflon  model  material,  lb/ft 

I  .hr 

C  „„„  =  Specific  heat  of  Teflon  model  material,  BTU/lb-°R 

p,  iEr 
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1 


=  Thermal  conductivity  of  Teflon  model  material,  BUT/ft-hr-uR 

TEr 

t  =  Time  for  initial  exposure  to  the  flow,  hr 
A  constant  value  of  2. 190  BUT/ft2  -hr*  5-°R  was  used  for^P  TEp  Cp^  TEp 

and  t  =  0  was  established  at ,  85  sec  before  the  model  reached  tunnel  centerline. 
Thermocouples  embedded  in  the  Teflon  flat  plate  insert  enabled  measurement  of  model 
material  temperature  prior  to  injection.  The  Teflon  booster  and  orbiter  models  were 
not  so  equipped  and  initial  material  temperature  was  assumed  as  ambient.  Accuracy 
of  the  phase  change  paint  technique  is  discussed  in  Ref.  7. 

c.  Pressure  Data 


The  model  surface  pressures  were  measured  by  two  types  of  transducers.  All 
odd-numbered  ports  were  connected  to  15  paid  transducers  while  the  even-numbered 
ports  were  connected  to  paralleled  1  psid  and  15  psid  transducers.  For  the  even- 
numbered  ports,  all  pressures  below  1  psia  were  read  on  the  1  psid  transducers  to 
give  accurate  results  for  the  undisturbed  pressures.  Uncertainties  are  given  below: 


Transducer 


Uncertaint 


1  psid  ±  0. 001  psia 

15  psid  ±  0. 003  psia  or 

1/2%  of  reading,  whichever 
is  greater. 


d.  Data  Output 

The  test  output  consists  of  tabulations,  data  plots  and  photographs  for  each  run. 

All  runs  are  denoted  by  a  group  number.  The  total  data  output  is  contained  in  22  volumes 
issued  by  AEDC,  and  is  catalogued  in  detail  in  Ref.  6. 
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Section  in 

TEST  DATA  OBTAINED 

Pressure  and  heating  rate  distributions  and  profile  flow  photographs  were  obtained 
for  the  flat  plate  and  booster  models  with  various  shock  genera  jrs  for  the  test  conditions 
given  in  paragraph  II— 2.  As  described  in  the  previous  section,  two  methods  were  used  to 

1  i 

obtain  the  heating  rate  distributions:  thermocouples  on  thin-skin  models,  and  phase  change 

I 

paints  on  teflon  models.  Both  methods  use  initial  heating  of  the  model.  Thus,  the  heat 
transfer  data  obtained  are  for  "cold  wall"  conditions.  The  static  pressure  measurements, 
conversely,  require  several  minutes  to  stabilize.  During  this  stabilization  time,  the  model 

i 

heats  and  the  wall  temperature  approaches  the  adiabatic  wall  temperature.  Therefore,  the 
pressure  data  are  for  essentially  adiabatic  wall  conditions,  whereas  the  heat  .transfer  data 

i 

are  for  cold  wall  conditions.  Of  course,  this  difference  in  conditions  is  undesirable  because 
wall  temperature  affects  both  boundary  layer  transition  and  also  separation.  It  was  expe¬ 
dient,  however,  to  fabricate  the  models  for  standard  data  measuring  techniques  and  thereby 
avoid  the  complexities  of  either  cooling  the  pressure  models  or  using  steady  state  heat  trans¬ 
fer  gauges. 

! 

The  data  obtained  on  the  flat  plate  and  booster  models  with  the  various  generators  are  • 
indicated  in  this  section.  The  sample  runs  chosen  ure  used  to  describe  features  of  different 
types  of  interaction  flows  and  the  effects  on  the  resulting  pressure  and  heat  transfer  distri¬ 
butions.  Finally,  a  complete  data  log  presents  a  summary  of  all  data  obtained. 

1.  SHARP  PLATE 

The  shock  generators  used  and  the  types  of  data  obtained  on  the  sharp  leading  edge  flat 
plate  are  indicated  in  Table  II.  The  tunnel  stagnation  pressure  (pQ)  and  nominal  location  of 
the  generated  shock  on  the  plate  surface  (x.)  are  listed  in  the  first  two  columns1.  The  remain¬ 
ing  columns  indicate  the  shock  generators  and  the  type  of  data  obtained  (p  =  pressure,  T  = 
thermocouple  heating  rates,  P  =  temperature-sensitive  paint  heating  rates).  The  row  with 
"trip"  in  the  first  column  indicates  the  data  obtained  with  the  boundary  layer  trip  installed  on 
the  flat  plate. 

Pressure  distributions  were  obtained  for  nearly  all  cases  where  thermocouple  heating 
rate  distributions  were  obtained.  However,  only  limited  temperature  sensitive  paint  data 
were  obtained.  Profile  flow  photographs  were  taken  for  all  test  configurations;  shadowgraphs 
on  thermocouple  runs  and  schlierens  on  pressure  runs. 
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a.  Undisturbed  Flow 

Undisturbed  (no  shock  generator)  pressure  and  heating  rate  distributions  on  the  flat 
plate  surface  were  measured  for  all  tunnel  stagnation  pressure  levels.  These  were  used  as 

I 

reference  conditions  to  determine  the  pressure  rises  and  heating  amplifications  caused  by 
the  generated  shock  wave  boundary  layer  interactions. 

The  heat  transfer  distributions,  obtained  on  the  thin-wall  thermocouple  model  with  no 
houndary  layer  trip,  exhibited  an  anomalous  behavior  near  the  trailing  edge  of  the  plate.  This 
anomalous  behavior,  which  was  most  pronounced  at  a  tunnel  flow  stagnation  pressure  of  400 
psia,  is  evidenced  by  the  heat  transfer  data  shown  in  Fig.  22.  Although  there  is  considerable 
scatter,  'particularly  for  the  lower  heating  rates,  the  distribution  dips  (at  x  =  23. 5  in.)  and 
then  rises  at  the  trailing  edge.  We  attribute  this  to  separation  and  reattachment  of  the  bound¬ 
ary  layer  due  to  a. slight  warping  of  the  thin  wall  surface  near  the  trailing  edge  of  the  plate. 

The  slight  warping  can  also  be  ascertained  by  examining  carefully  the  shadowgraph 

» 

photographs  taken  during  the  same  tunnel  run  (Fig.  23).  The  triangular  shapes  on  the  lower 
surface  of  the  model  are  photo  reference  guides  located  8. 0  and  22. 0  in.  downstream  of  the 
leading  edge.  The  square  grid  on  the  tunnel  windows  has  a  spacing  of  4. 5  in.  In  addition  to 
showing  the  slight  concave,  curvature  of  the  plate  trailing  edge,  these  shadowgraphs  clearly 
show  a  displaced  white  line  above  the  plate  surface,  indicating  that  the  boundary  layer  was 
laminar  over  the  entire  extent  of  the  plate.  ^ 

'The  measured  heating  rates  for  tunnel  flow  stagnation  pressures  of  850  psia  and  200 
psia  also  exhibit  an  anomalous  dip  near  the  plate  trailing  edge  (Fig.  24),  but  not  as  severe 
as  for  the  400  psia  case  mentioned  above.  When  the  sharp  plate  receiver  was  tested  with 
shock  generators,  it  was  found  that  the  great  majority  of  peak  interaction  heating  rates  were 
measured  .well  upstream  of  the  regions  that  appeared  separated  on  the  undisturbed  plate,  and 
therefore  would  be  unaffected  by  the  separation.  There  were  five  cases  where  the  peak  in¬ 
teraction  heating  rates  were  measured  upstream  of,  but  in  proximity  to,  the  beginning  of  the 
separated  regions.  It  is  believed  that  in  these  cases  as  well,  the  measured  peak  heating 
rates  yvere  unaffected  by  warping  of  the  plate  trailing  edge,  primarily  because  the  peaks 
occurred  upstream  of  the  separated  regions,  but  also  because  these  data  were  not  inconsis¬ 
tent  with  all  other  sharp  plate  data.  For  these  reasons,  the  smoothly  faired  curves  of 

i  I 

^As  described  by  Chapman,  Kuehn  and  Larson  (Ref.  8),  a  white  line  appears  above  laminar 
boundary  layers  in  shadowgraph  photographs.  This  line  remains  parallel  to  the  plate  surface 
while  the  boundary  remains  laminar,  converges  to  the  surface  in  the  boundary  layer  trans¬ 
ition  region,  and  disappears  when  the  flow  is  fully  turbulent. 

I 
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Figs.  22  and  24,  were  used  to  define  reference  levels  of  undisturbed  heating.  There  was 
one  case  (Group  60;  1/2  in.  dia.  sphere;  pQ  =  400  psia;  x.  =  22  in.)  where  the  peak  heating 
rate  was  measured  at  a  point  within  the  region  affected  by  the  warping  of  the  plate,  and  this  is 
noted  in  the  data  summary  of  Table  HI.  It  has  not  been  determined  if  the  peak  heating  rate 
in  this  case  was  affected  by  the  warped  trailing  edge. 

Although  a  Mach  wave  emanating  from  the  nose-plate  junction  is  visible  in  the  shadow¬ 
graph  photographs  (joint  at  x  =  2. 5  in. ,  see  Fig.  23),  this  joint  was  quite  smooth  and  did  not 
affect  either  the  heat  transfer  or  pressure  distributions.  The  pressure  model  remained 
quite  smooth  (no  local  warping  of  the  surface),  and  there  were  no  anomalies  in  the  undis¬ 
turbed  pressure  distributions  (see  Fig.  25). 


Schlieren  and  shadowgraph  flow  photographs  (Fig.  26)  show  no  boundary  layer  trans¬ 
ition  on  the  plate  surface  even  for  the  highest  tunnel  flow  stagnation  pressure  (850  psia).  In 
order  to  obtain  turbulent  boundary  layer  interaction  data  on  the  sharp  flat  plate,  it  was  nec¬ 
essary  to  use  boundary  layer  trips  (described  in  the  previous  section).  As  indicated  in 
Table  II,  the  trips  were  used  on  the  sharp  flat  plate  only  for  the  highest  tunnel  flow  stagna¬ 
tion  pressure  (850  psia). 


Schlieren  and  shadowgraph  photographs  of  the  tripped  boundary  layer  (taken  during  the 
pressure  and  heat  transfer  tunnel  runs,  respectively)  are  shown  in  Fig.  27.  The  corres¬ 
ponding  pressure  and  heat  transfer  rate  distributions  are  shown  in  Fig.  28a  and  b,  respec¬ 
tively.  The  curves  faired  through  the  data  were  used  as  undisturbed  reference  conditions 
for  the  tripped  boundary  layer  interactions.  Heating  rate  distributions  for  laminar  and  tur¬ 
bulent  boundaries  on  a  sharp  flat  plate  with  constant  wall  temperature  were  calculated  using 
Eckert's  reference  temperature  method  (Ref  9): 


hLAM 


for  laminar  boundary  layers,  and 


k+  — 

bruRB  =  0.0296  ”(Pr*)  3  (Re^)  °*8 


for  turbulent  boundary  layers.  The  reference  temperature  is: 

T*  =  4  <T  +  T  )  +  0.22  (T  -  T  ) 

*  2  v  w  e '  '  r  e’ 
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where 

T^  =  wall  temperature 

Tg  =  temperature  at  boundary  layer  edge 

Tf  =  recovery  temperature 

h  =  heat  transfer  coefficient 

Pr  =  Prandtl  number 

Re  =  Reynolds  number  based  on  x 

A* 

k  =  gas  thermal  conductivity 

and  *  conditions  are  based  on  the  reference  temperature.  Heating  rates  predicted  using  this 
method  are  also  shown  in  Fig.  28b.  Reference  conditions  for  the  heat  transfer  data  are 
listed  in  Table  IV. 

b.  Laminar  Upstream  of  Interaction 

In  all  cases  without  the  trip,  the  boundary  layer  was  initially  laminar  upstream  of  the 
interaction  caused  by  the  generator  shocks.  The  generator  shock  waves  were  usually  strong 
enough  to  separate  the  boundary  layer,  and  frequently  were  sufficiently  strong  to  cause 
boundary  layer  transition  prior  to  reattachment  (Ref.  10). 

Heating  rate  and  pressure  distributions  for  a  shock  wave  generated  by  a  10°  wedge 
incident  on  the  flat  plate  surface  at  a  nominal  value  of  x  =  22  in. ,  are  plotted  in  Fig.  29. 
These  distributions,  obtained  for  a  tunnel  stagnation  pressure  of  400  psia,  are  compared 
with  the  undisturbed  distributions  for  the  same  tunnel  pressure.  The  heat  transfer  starts  to 
dip  below  the  undisturbed  value  approximately  11-1/4  in.  downstre;  of  the  plate  leading 
edge.  The  schlieren  photograph  in  Fig.  30  confirms  this  as  the  sta  ;  of  a  separated  flow 
region.  Unfortunately,  the  shadowgraph  does  not  show  the  start  of  separation  because  the 
model  was  I  ,cated  further  downstream  and  the  separation  location  was  between  the  two  tun¬ 
nel  viewing  windows. 

In  this  case,  the  shock  wave  emanating  from  the  separation  location  impinges  on  the 
wedge,  is  reflected,  and  merges  with  the  wedge  shock  (evident  in  both  the  shadowgraph  and 
schlieren  in  Fig.  30).  The  flow  in  the  impingement  region  is  quite  complex  in  these  cases. 
Nevertheless,  the  peak  heating  rate  and  pressure  values  occur  approximately  at  the  chosen 
nominal  impingement  location,  22  in.  downstream  of  the  leading  edge.  The  peak  heating 
rate  is  more  than  60  times  the  undisturbed  value;  the  peak  pressure  is  nearly  30  times  as 
large  as  the  undisturbed  pressure. 

The  pressure  measurements  were  made  by  reading  16  pressure  transducers  at  8 
different  valve  positions.  Since  several  minutes  were  required  for  the  pressure  readings  at 
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any  one  valve  position  to  stabilize,  only  the  portion  of  the  pressure  distributions  in  the 
region  of  impingement  was  recorded  during  most  tunnel  runs.  It  was  therefore  possible  to 
obtain  more  peak  pressure  rises  than  would  have  been  possible  if  complete  pressure  distri¬ 
butions  had  been  taken.  This  was  desirable  in  order  to  obtain  a  maximum  number  of  data 
points  for  the  correlations  of  peak  pressure  and  heating  described  in  the  following  section. 
Of  course,  incomplete  pressure  distributions  make  it  more  difficult  to  correlate  other 
parameters  such  as  separation  distance  and  wall  temperature  effects. 

Finally,  also  evident  in  Fig.  30,  is  a  lip  shock  emanating  from  the  trailing  edge  of 
the  wedge,  immediately  downstream  of  the  expansion  fan.  As  described  by  Hama  (Ref.  11), 
the  flow  around  the  expansion  corner  (the  wedge  trailing  edge  in  our  case)  overexpands 
appreciably  and  then  is  recompressed  through  the  lip  shock  and  separates  from  the  wedge 
base.  In  certain  instances,  this  lip  shock  impinged  on  the  plate  surface  and  led  to  locally 
increased  pressures  and  heating  rates  on  the  plate  surface. 

The  shock  wave  emanating  from  the  boundary  layer  separation  location  on  the  plate 
surface  was  sufficiently  strong  in  some  cases  to  cause  boundary  layer  separation  from  the 
wedge  surface.  This  was  the  case,  for  example,  for  the  15  deg  wedge  shock  impinging  at 
x.  =  5  in.  for  a  tunnel  flow  stagnation  pressure  of  75  psia.  The  heat  transfer  and  pressure 
distributions  (Fig.  31),  and  the  corresponding  shadowgraph  and  schlieren  flow  photographs 
(Fig.  32),  indicate  laminar  boundary  layer  separation  essentially  from  the  leading  edge  of 
the  plate  with  an  accompanying  separation  shock.  This  shock  merged  with  the  plate  bow 
wave  and  impinged  near  the  ti-ailing  edge  of  the  wedge,  causing  the  laminar  boundary  layer 
on  the  wedge  to  separate.  For  the  "cold  wall"  case  (shadowgraph),  separation  occurred 
downstream  of  the  wedge  mid-chord  (cold  walls  general! y  delay  separation);  for  the  "hot 
wall"  case  (schlieren),  separation  occurred  upstream  of  the  wedge  mid-chord.  In  both 
cases,  the  wedge  separation  shock  merged  with  the  shock  wave  generated  by  the  wedge 
before  interacting  with  the  plate  boundary  layer.  The  peak  pressure  was  measured  on  the 
plate  surface  at  x«s*4. 8  in. ,  and  the  peak  heating  rate  on  the  plate  was  measured  at 
x  5.2  in. 

In  both  cases,  the  plate  shock  is  reflected  from  near  the  wedge  trailing  edge  and 
results  in  additional  local  peaks  (most  notably  in  the  pressure  data)  in  the  heating  rate 
ajid  pressure  distributions  on  the  plate  surface. 

When  the  wedge  was  moved  further  aft,  on  occasion  the  plate  bow  wave  impinged 
essentially  at  the  wedge  leading  edge  and,  particularly  for  the  small  wedge  angle  (1. 5  deg), 
considerably  strengthened  the  shock  wave  generated  by  the  wedge  (Figs.  33  and  34).  For 
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the  1. 5  deg  wedge  with  a  nominal  x.  =  11  in. ,  and  a  tunnel  flow  stagnation  pressure  of  850 
psia,  the  wedge  generated  shock  wave  impinged  on  the  plate  boundary  layer,  was  reflected, 
impinged  on  the  wedge  near  the  trailing  edge,  was  again  reflected  and  again  impinged  on  the 
plate  boundary  layer.  This  caused  the  multiple  peaks  in  the  heat  transfer  and  pressure  dis¬ 
tributions  evident  in  Fig.  33. 
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Heat  transfer  and  pressure  distributions  along  the  plate  centerline,  and  the  corres¬ 
ponding  shadowgraph  and  schlieren  flow  photographs,  for  a  10  deg  wedge  shock  impinging  at 
a  nominal  location  22  in.  downstream  of  the  leading  edge,  are  shown  in  Fig.  35  and  36  for 
a  tunnel  stagnation  pressure  of  850  psia.  Data  distributions  and  flow  photographs  for  a  15 
deg  wedge  shock,  for  the  same  tunnel  flow  pressure  and  nominal  impingement  location,  are 
shown  in  Fig.  37  and  38.  The  shadowgraph  photographs  for  both  wedges  indicate  a  separa¬ 
tion  location  further  downstream  than  that  shown  in  the  corresponding  schlieren  photographs. 
The  "cold  wall"  delayed  separation  of  the  boundary  layer  from  the  plate  surface.  This  effect 
is  also  evidenced  by  the  centerline  data  distributions  in  Fig.  37.  The  heat  transfer  data  first 
diverge  from  the  undisturbed  values  approximately  10. 5  in.  downstream  of  the  leading  edge. 
Whereas  the  pressure  data,  obtained  on  a  "hot  wall"  model,  start  to  increase  above  the  un¬ 
disturbed  pressure  values  approximately  9. 5  in.  downstream  of  the  leading  edge. 

The  pressure  rises  to  a  laminar  separation  plateau  level  that  is  well  approximated  by 
Hill's  (Ref  12)  correlation: 
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where  pp^  is  the  laminar  plateau  pressure  level,  subscript  1  refers  to  undisturbed  con¬ 
ditions,  and  Re  is  the  Reynolds  number  based  on  undisturbed  conditions  and  distance  from 
sep 

the  leading  edge  to  the  separation  location.  The  value  obtained  from  Eq  5,  ppy  =  1.7  p^, 
as  shown  in  Fig.  37,  is  just  slightly  below  the  measured  values. 


The  measured  pressure  amplifications  were  compared  with  calculated  amplifications 
and  a  number  of  significant  differences  and  trends  were  noted.  These  are  discussed  below. 


(1)  Calculated  Pressure  Ratios 

For  wedge  and  cone  generators,  oblique  shock  relations  (Ref.  13)  were  used  to  obtain 
the  shock  angle  and  the  inviscid  properties  behind  the  shock.  At  this  point,  a  first  order 
correction  for  the  effect  of  the  generator  boundary  layer  disp'  .cement  thickness  was  made. 
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Using  the  computed  downstream  properties  and  simplified  boundary  layer  equations  (Ref  14), 
the  generator  boundary  layer  displacement  thickness  was  computed,  and  the  corrected  cone 
or  wedge  angle  was  determined  (Fig.  39). 

Properties  downstream  of  the  generated  shock  were  recomputed  using  the  corrected 
generator  angle.  It  was  assumed  that  along  the  plate  center  line  it  would  be  acceptable  to 
treat  the  reflected  shock  as  a  planar  shock  for  the  case  of  a  cone  generator,  as  well  as  for 
the  wedge  generators.  Oblique  shock  relations  were  then  used  to  determine  the  properties 
behind  the  reflected  shock. 

For  the  hemisphere-cylinder  generators,  it  was  first  necessary  to  define  the  shape  of 
the  incident  shock.  For  this  purpose  Lukasiewicz  (Ref  15)  blast  analogy  for  hypersonic  flow 
was  used.  The  local  shock  angle  at  the  centerline  impingement  point  was  computed,  and  the 
wave  at  this  location  was  treated  as  a  section  of  an  oblique  shock.  The  remainder  of  the 
computation  then  becomes  identical  to  the  wedge  calculation,  although  no  boundary  layer 
correction  was  required. 

(2)  Measured  Pressure  Ratios  -  Wedge  Generators 

The  flow  fields  observed  during  testing  were  significantly  more  complex  than  the 
idealized  flow  field  shown  in  Fig.  39,  resulting  in  differences  between  measured  and  pre¬ 
dicted  pressure  amplifications.  The  general  trends  observed  in  the  wedge/sharp  flat  plate 
data  are  indicated  in  Fig.  40,  where  measured  pressure  ratios  are  compared  with  calculated 
ratios  for  testing  at  850  psia  and  75  psia.  These  data  are  presented  at  "r.  minal  impinge¬ 
ment  locations",  which  are  the  idealized  locations  calculated  for  the  simplified  flowfield. 
Although  these  locations  are  referred  to  throughout  the  report  for  simplicity,  the  actual 
locations  for  impingement  with  the  plate,  as  well  as  peak  amplification,  differ  from  the 
nominal  values,  as  seen  in  Table  HI. 

As  described  previously,  for  the  11-  and  22-inch  impingement  locations,  there 
are  additional  shock  waves  that  impinge  on  the  flat  plate  and  increase  the  pressure  ratio 
above  that  value  predicted  for  a  simple  shock  generated  by  the  wedge.  Four  phenomena, 
each  leading  to  an  additional  shock  wave,  are:  1)  reflection  of  the  plate  leading  edge 
shock  by  the  wedge  surface,  2)  reflection  of  a  shock  wave  emanating  from  the  start  of  a 
separated  flow  region  on  the  plate  by  the  wedge  surface,  3)  a  shock  wave  emanating 
from  a  separated  flow  region  on  the  wedge  surface,  and  4)  a  wedge  lip  shock.  In  each  case, 
the  additional  shock  reinforces  the  simple  wedge  generated  shock  wave  and  results  in  mea¬ 
sured  pressure  amplifications  on  the  plate  surface  that  are  substantially  larger  than  those 
predicted  for  the  simple  wedge  shock  (Fig.  39). 
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However,  at  the  forward  (5  inch)  impingement  location,  the  measured  pressure  am¬ 
plifications  agree  relatively  well  with  the  calculated  amplifications.  To  obtain  a  forward 
impingement,  the  generator  was  positioned  with  its  leading  edge  upstream  of  the  plate 
leading  edge  shock.  With  the  generator  in  this  position,  the  forward  portion  of  the  plate 
tended  to  be  free  of  extraneous  shocks  (Figure  41),  with  separation  shocks  a  common  ex¬ 
ception  (Figure  42).  As  seen  in  these  figures,  the  reflection  of  the  plate  leading  edge  shock 
from  the  wedge  strikes  the  plate  well  downstream  of  the  nominal  impingement  point.  The 
separation  shocks  were  either  not  reflected  (Figure  42),  or  if  reflected,  they  also  struck 
the  plate  well  aft  of  the  nominal  impingement  point.  The  result  is  that  with  the  wedge  gener¬ 
ator  in  the  forward  position,  there  was  no  impingement  of  extraneous  shocks  in  the  interac¬ 
tion  region  or  upstream  of  this  region.  The  net  pressure  amplification  (p^g^/py^)  there¬ 
fore  remained  approximately  unchanged  from  the  value  predicted  for  the  idealized  flow  field. 

(3)  Three-Dimensional  Generators-Cones  and  Hemisphere  Cylinders 

Measured  pressure  amplifications  for  the  cone  and  hemisphere  cylinder  generators 
were  significantly  lower  than  predicted  (Fig.  43),  and  therefore  also  below  the  wedge  gener¬ 
ator  data.  It  is  believed  that  this  trend  is  a  result  of  the  highly  three-dimensional  nature  of 
the  flow  field  downstream  of  the  3D  generator  shock. 

The  strength  of  the  reflected  3D  generator  shock  continuously  decreases  away  from  the 
plate  centerline.  The  resulting  divergence  of  the  flow  downstream  of  the  reflected  shock 
provides  relief  which  reduces  the  pressure  amplification.  The  divergent  nature  of  the  flow 
downstream  of  the  3D  generator  shock  is  illustrated  in  Figure  44  which  compares  the  lateral 
heat  transfer  and  pressure  distributions  for  wedge  and  cone  generators  of  approximate¬ 
ly  equal  calculated  shock  strengths.  It  should  be  noted  that  the  transverse  data  were  ob¬ 
tained  downstream  of  the  point  of  peak  amplification. 

As  a  result  of  three  dimensional  flow  relief,  the  3D  generators  did  not  produce  the 
large  separated  regions  in  the  plate  boundary  layer,  and  the  resulting  separation  shocks 
that  were  observed  with  the  wedge  generators.  This  is  seen  in  Figures  45  and  46  where  the 
interactions  caused  by  a  wedge,  cone,  and  hemisphere-cylinder  of  approximately  equal  shock 
strengths  are  compared.  The  wedge  generator  shock  produced  an  extensive  separated  re¬ 
gion  (Fig.  46a).  The  accompanying  separation  shock  is  reflected  by  the  generator  and  merges 
with  the  generator  shock  prior  to  impinging  on  the  plate  boundary  layer  (Fig.  45a).  There  is 
no  visible  evidence  in  the  shadowgraphs  of  Figures  45b  and  45c  that  the  cone  and  hemisphere 
cylinder  generators  caused  separation  of  the  plate  boundary  layer.  The  heat  transfer  dis¬ 
tributions  for  the  3D  generators  (Figs  46b  and  46c),  however,  indicate  possible  small  regions 
of  separation. 
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The  three  dimensional  nature  of  the  interaction  heating  distribution  caused  by  the 
2. 0  in.  diameter  hemisphere-cylinder  is  further  evidenced  by  photographs  taken  during  the 
temperature-sensitive  paint  runs,  as  shown  in  Fig.  47.  The  particular  photograph  shown 
in  this  figure  was  taken  26. 92  secs  after  the  model  first  entered  the  tunnel  flow.  The 
paint  chosen  for  this  run  had  a  phase  change  temperature  of  710°R,  which  was  the  tempera¬ 
ture  along  the  boundary  of  the  dark  crescent  shaped  region  on  the  plate  surface  at  the  time 
the  photograph  was  taken.  The  corresponding  heating  rate  at  the  boundary  (hT  =  0. 21 
agrees  with  the  thermocouple  data  shown  in  the  figure. 

c.  Turbulent  Upstream  of  Interaction 

As  noted  earlier  in  this  section,  even  at  the  highest  tunnel  flow  stagnation  pressure 
(850  psia)  it  was  necessary  to  use  boundary  layer  trips  in  order  to  obtain  turbulent  boundary 
layers  on  the  sharp  leading  edge  flat  plate.  In  order  to  have  as  high  a  Reynolds  number  as 
possible,  the  nominal  impingement  location  was  chosen  near  the  trailing  edge  of  the  plate  for 
all  tripped  boundary  layer  cases  (see  Table  II).  A  comparison  of  the  undisturbed  flows  over 
the  sharp  leading  edge  plate,  without  and  with  the  boundary  layer  trip,  is  shown  by  the  flow 
photographs  in  Figs.  26  and  27.  The  shadowgraph  photographs  were  obtained  on  a  colder 
wall  than  the  schlieren  photographs,  which  is  expected  to  delay  transition.  The  heat  trans¬ 
fer  data,  shown  in  Fig.  28,  indicate  that  the  tripped  boundary  layer  became  hilly  turbulent 
(end  of  transition)  approximately  18  in.  downstream  of  the  leading  edge.  For  the  pressure 
data,  obtained  on  hot  wall  models,  transition  is  expected  to  have  occurred  somewhat  further 
upstream. 

Some  effects  of  the  boundary  layer  trip  on  the  interaction  caused  by  a  10  deg  wedge 
shock  generator  can  be  ascertained  by  comparing  the  untripped  data  of  Fig.  35  and  36  with 
the  tripped  boundary  layer  data  shown  in  Figs.  48  and  49.  For  the  colder  wall  (heat  transfer 
data),  Fig.  35a  indicates  separation  approximately  11  in.  downstream  of  the  leading  edge, 
whereas  for  the  tripped  boundary  layer,  Fig.  48a  indicates  separation  approximately  20  in. 
downstream  of  the  leading  edge.  Similarly,  the  pressure  distributions  also  indicate  that 
separation  is  delayed  by  the  trip.  Thus,  as  expected,  separation  is  much  less  extensive  for 
the  tripped  turbulent  boundary  layer  than  for  the  initially  laminar  boundary  layer  for  both 
the  "cold"  and  "hot  wall"  models. 

The  tripped  boundary  layers,  with  relatively  small  separated  regions,  did  not  exhibit 
the  strong  separation  shocks  which  were  characteristic  of  the  interactions  between  wedge 
generator  shocks  and  laminar  boundary  layers.  As  a  result,  the  measured  pressure  ampli¬ 
fications  for  wedge  generators  and  tripped  boundary  layers  were  lower  than  the 
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amplifications  measured  for  the  untripped  boundary  layers,  and  in  general  showed  better 
agreement  with  predicted  values  (Fig.  50). 

The  15  deg  wedge  shock  generator  was  used  at  two  different  heights  above  the  flat 
plate  for  the  tripped  boundary  layer.  Heat  transfer  and  pressure  distributions  for  the  lower 
height,  shown  in  Fig.  51,  indicate  separation  approximately  19-1/2  in.  downstream  of  the 
leading  edge  for  the  cold  wall  (heating  data),  whereas  separation  occurred  less  than  15  in. 
downstream  of  the  leading  edge  for  the  hot  wall  (pressure  data).  The  pressure  distribution 
exhibits  a  plateau  region  indicating  that  the  boundary  layer  was  initially  laminar  at  the  be¬ 
ginning  of  the  separated  flow  region.  The  schlieren  photographs  in  Fig.  52  support  the  con¬ 
clusion  that  the  extent  of  separation  was  considerably  larger  for  the  pressure  run  than  for 
the  heat  transfer  run.  At  the  greater  generator  height,  the  data  distributions  and  flow 
photographs  indicate  turbulent  separation  for  both  the  heat  transfer  and  the  pressure  runs 
(Figs.  53  and  54). 

Healing  rate  distributions  on  the  sharp  leading  edge  flat  plate  with  the  boundary  layer 
trip  were  obtained  with  both  thermocouples  and  temperature-sensitive  paint.  Thermocouple 
heating  data  and  the  pressure  distribution  resulting  from  the  interaction  caused  by  a  7-1/2 
deg  cone  are  presented  in  Fig.  55.  They  indicate  a  relatively  small  region  of  increased 
heating  and  pressure  starting  approximately  20  in.  downstream  of  the  leading  edge.  Flow 
photographs  for  this  configuration  are  shown  in  Fig.  56.  Photographs  showing  temperature- 
sensitive  paint  results,  and  the  corresponding  heat  transfer  rates,  are  presented  in  Fig.  57. 
Times  are  measured  from  when  the  model  first  enters  the  tunnel  flow.  For  this  run,  1. 23 
seconds  were  required  from  time  zero  until  the  model  was  on  the  tunnel  centerline.  The 
phase  change  temperature  of  the  paint  chosen  was  710°R.  In  each  frame,  this  is  the  tem¬ 
perature  along  the  boundary  of  the  melted  (dark  colored)  and  unmelted  (light  colored)  paint 
regions.  This  boundary  also  marks  a  line  of  constant  heating  rate;  the  heating  rates  are 
lower  in  the  light  colored  region  and  higher  in  the  dark  colored,  crescent  shaped,  region  on 
the  plate  surface. 

Because  a  spherical  shock  wave  is  attenuated  by  expansion  waves  and  becomes  weaker 
away  from  the  centerline,  one  might  expect  a  conical  shock  wave  to  cause  a  greater  dis¬ 
turbance  outboard  of  the  centerline  than  that  caused  by  the  spherical  shock  wave,  for  equal 
shock  strengths  on  the  centerline  and  approximately  equal  peak  heating  rates.  The  photo¬ 
graphs  shown  in  Fig.  58a  bear  this  out  For  equal  heating  rates,  the  conical  shock  inter¬ 
action  with  the  turbulent  boundary  layer  is  more  extensive  than  the  spherical  shock  inter¬ 
action.  However,  for  laminar  boundary  layers  (Fig.  58b),  the  conical  and  spherical  shock 
interactions  extend  equal  distances  outboard  of  the  centerline. 
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2.  BLUNT  PLATE 


Interaction  data  on  the  blunt  leading  edge  flat  plate  were  obtained  at  three  tunnel  flow 
stagnation  pressure  levels  (850,  400  and  75  psia)  for  two  nominal  impingement  locations 
(x.  =  7  in.  and  22  in, ).  The  particular  shock  generators  used  are  indicated  in  Table  V.  As 
indicated  in  the  table,  no  temperature  sensitive  paint  data  were  obtained  for  the  blunt 
leading  edge  flat  plate. 

As  with  testing  of  the  sharp  leading  edge  plate,  shadowgraph  flow  photographs  were 
obtained  during  the  initial  heating  of  the  model  on  thermocouple  runs  and  schlieren  photo¬ 
graphs  were  obtained  during  the  pressure  data  runs  when  the  model  wall  was  hot. 

a.  Undisturbed 

The  blunt  leading  edge  led  to  a  substantial  favorable  pressure  gradient  over  the  entire 
flat  plate  surface  as  illustrated  by  Fig.  59,  for  a  tunnel  flow  stagnation  pressure  of  850  psia. 
The  corresponding  undisturbed  heating  distribution  indicates  chat  boundary  layer  transition 
started  approximately  14  in.  downstream  of  the  blunt  leading  edge  for  850  psia.  The  shadow¬ 
graph  and  schlieren  flow  photographs  for  these  tunnel  runs  are  shown  in  Fig.  60.  Comparing 
the  flow  photographs  and  heat  transfer  data,  for  the  highest  tunnel  flow  pressure,  it  appears 
that  the  boundary  layer  was  transitional  over  the  downstream  portion  of  the  flat  plate. 

Undisturbed  distributions  at  the  lower  tunnel  pressures  indicated  entirely  laminar 
boundary  layer  flow  over  the  flat  plate  surface.  Curves  were  faired  through  all  the  un¬ 
disturbed  data  distributions,  see  Fig.  61,  and  used  as  reference  conditions  for  the  down¬ 
stream  (Xj«22  in.)  and  upstream  (x.*»7  in.)  interactions. 

b.  Downstream  Interactions 

Centerline  heating  and  pressure  distributions,  for  a  10  deg  wedge  generating  a  shock 
impinging  at  a  nominal  location  approximately  22  in.  downstream  of  the  leading  edge,  are 
shown  in  Fig.  62  for  a  tunnel  flow  stagnation  pressure  of  850  psia.  The  corresponding 
shadowgraph  and  schlieren  photographs  are  in  Fig.  63.  The  cold  wall  heating  rate  dis¬ 
tribution  initially  dips  below  the  undisturbed  distribution,  approximately  16  1/2  in.  down¬ 
stream  of  the  leading  edge,  and  then  attains  a  peak  heating  amplification  at  approximately 
21  1/2  in.  The  pressure  distribution  rises  upstream  of  15  in.  to  a  plateau  level  and  then 
attains  a  peak  value  21  in.  downstream  of  the  leading  edge. 

Similar  distributions  and  flow  photographs,  but  for  the  lowest  tunnel  flow  pressure 
(75  psia),  are  given  in  Fig.  64  and  65.  Separation  is  considerably  more  extensive  for  the 
low  tunnel  pressure  than  for  the  high  one.  The  heat  transfer  rates  drop  below  the  undisturbed 
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values  11  in.  downstream  of  the  leading  edge,  and  attain  a  oeak  value  at  19  in.  The  limited 
pressure  distribution  indicates  separation  well  upstream  of  15  in.  and  a  peak  value  just 
downstream  of  19  in.  A  second,  lower,  peak  is  apparent  in  both  the  heating  and  pressure 
distributions.  This  is  probably  caused  by  additional  shock  waves  being  reflected  by  the 
wedge  surface  and  impinging  on  the  plate.  However,  the  density  at  these  low  tunnel  pres¬ 
sures  is  insufficient  to  discern  reflected  shocks  in  the  flow  photographs. 

c.  Upstream  Interactions 

Heating  and  pressure  distributions  for  the  upstream  impingement  location  are  shown 
in  Fig.  66  for  a  10  deg  wedge  for  850  psia,  the  corresponding  flow  photographs  are  shown 
in  Fig.  67.  These  photographs  clearly  show  the  strong  and  complex  interaction  flow 
between  the  wedge  and  plate.  The  heating  rate  is  an  order  of  magnitude  higher  than  the 
undisturbed  value  at  the  start  of  the  instrumentation*  The  pressure  distribution  also  greatly 
exceeds  the  undisturbed  one;  even  at  the  most  forward  tap  location  (x  =  1  in.),  the  pressure 
is  four  times  larger  than  the  undisturbed  pressure.  The  pressure  level  remains  very  high 
until  near  the  wedge  trailing  edge,  and  then  drops,  indicating  an  accelerating  flow  downstream 
of  the  wedge  trailing  edge. 

Corresponding  data,  obtained  at  the  low  tunnel  pressure,  are  shown  in  Figs.  68  and 
69.  In  this  case  the  wedge  was  moved  aft  0. 6  in. ,  which  placed  the  wedge  leading  edge  just 
downstream  of  the  bow  shock  from  the  blunt  leading  edge  of  the  plate.  Again,  the  flow 
between  the  wedge  and  plate  is  quite  complex;  normal  shock  waves  are  visible  between  the 
wedge  and  plate  surfaces. 

When  the  wedge  was  pitched  to  15°  deg,  for  the  pressure  run,  the  model  blocked  the 
tunnel  flow.  It  was  necessary  to  increase  the  tunnel  flow  stagnation  pressure  to  100  psia, 
inject  the  model,  and  then  reduce  the  pressure  to  75  psia  in  order-  to  avoid  blocking  the 
tunnel  flow  for  the  pressure  run.  For  the  heat  transfer  run  for  this  model  configure  /  :n,  the 
tunnel  back  pressure  was  slightly  lower  and  the  tunnel  flow  started  during  a  second  injection 
of  the  model.  The  heating  and  pressure  distributions  for  the  15  deg  wedge,  and  the  corre¬ 
sponding  flow  photographs,  are  in  Fig.  70  and  71,  respectively.  The  heating  rate  peaks 
just  downstream  of  x  =  4  in.  to  a  value  14  times  higher  than  the  undisturbed  heating  rate. 

The  pressure  distribution  rises  at  x  =  4  in. ,  and  then  continues  to  rise  slowly  until  a  peak 
value  almost  15  times  as  large  as  the  undisturbed  pressure  is  attained  approximately  7  1/2 
in.  downstream  of  the  leading  edge.  [The  shadowgraph  photographs  in  Figs.  67  and  71  show 
a  shock  interaction  pattern  near  the  wedge  leading  edge  very  similar  to  one  first  described 
by  Edney  (Ref.  16  Fig.  6-11).  As  noted  by  Edney,  this  type  of  interaction  results  in  an 
impingement  on  the  wedge  surface  that  results  in  particularly  severe  heating.] 


23 


As  a  final  example  of  the  blunt  leading  edge  plate  data.  Figs.  72  and  73  show  the 
heating  distribution  measured  along  the  plate  centerline  and  the  corresponding  shadowgraph 
photograph  for  the  orbiter  shock  generator  parallel  to  the  plate  for  a  tunnel  stagnation 
pressure  of  850  psia.  The  peak  heating,  at  x  -  7  in.  is  almost  six  times  the  undisturbed 
heating  rate.  The  measured  heating  then  drops  to  a  level  approximately  4  times  higher 
than  the  undisturbed  level,  and  maintains  this  amplification  until  the  trailing  edge  of  the 
orbiter. 

3.  BOOSTER  AND  BOOSTER/ORBITER  CONFIGURATIONS 

Table  VT  presents  an  overview  of  the  data  obtained  on  the  booster  and  booster-orbiter 
configurations.  Again,  p  indicates  pressure  data,  T  indicates  thermocouple  heat  transfer 
data,  and  P  indicates  temperature-sensitive  paint  heat  transfer  data.  As  with  testing  of 
the  flat  plate  models,  heat  transfer  rates  were  obtained  during  the  initial  heating  of  the 
model,  using  both  the  thin-wall  thermocouple  and  temperature  sensitive  paint  techniques, 
and  therefore  the  heating  rates  correspond  to  a  relatively  ’’cold  wall"  condition.  Pressure 
data  were  obtained  for  model  wall  temperatures  approaching  the  adiabatic  wall  temperature 

Undisturbed  heating  rate  and  pressure  distributions  were  obtained  by  testing  the 
booster  alone  at  0  deg  angle  of  attack  and  the  orbiter  alone  at  0,  5  and  10  deg  angles  of 
attack.  As  indicated  in  Table  VI,  various  shock  generators  were  used  with  the  booster 
model.  The  orbiter  shock  generator  model  was  instrumented.  Heat  transfer  and  pressure 
distributions  were  measured  on  the  lower  surface  of  the  orbiter  as  well  as  the  upper  sur¬ 
face  cf  the  booster.  In  this  section  data  are  presented  for:  the  undisturbed  flows  over  the 
booster  and  orbiter  alone,  the  effects  on  the  booster  caused  by  the  interaction  flows,  and 
the  effects  on  the  orbiter  caused  by  the  interaction  flows. 

a.  Undisturbed  Flows 

Heat  transfer  rates  measured  on  the  upper  surface  of  the  booster  with  no  shock  gen¬ 
erator,  at  tunnel  flow  stagnation  pressures  of  850  and  75  psia,  are  indicated  in  Figs.  74 
and  75.  Two  streamwise  distributions  and  one  spanwise  distribution  are  shown  in  each 
figure.  One  streamwise  distribution  is  in  the  booster  centerplane  (row  1)  and  the  other  is 
along  the  surface  in  a  plane  at  30°  with  respect  to  the  centerplane  (row  2).  The  spanwise 
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distribution  is  on  the  wing  surface  at  station  x  =  15.6  in.  (see  Fig.  17).  The  off  centerline 
distributions  (row  2)  exhibit  an  increased  level  downstream  of  x  =  10  inches.  This  increase 
can  be  attributed  to  turbulence  introduced  by  the  shock  emanating  from  the  booster  canards. 
As  would  be  expected,  there  is  considerably  more  scatter  in  the  heat  transfer  rates  mea¬ 
sured  at  the  lo\,  v  stagnation  pressure. 

Streamwise  pressure  distributions  for  the  same  tunnel  pressures  are  given  in  Fig.  76 
in  the  centerplane  (row  1)  and  in  a  plane  at  30°  with  respect  to  the  centerplane  (row  3).  The 
pressures  do  not  exhibit  a  sizable  jump  downstream  of  the  booster  canards;  the  distributions 
along  the  two  rows  are  essentially  identical  for  each  tunnel  flow  pressure.  However,  booster 
pressures  at  the  75  psia  tunnel  pressure  are  slightly  higher  than  those  at  the  850  psia  tunnel 
pressure,  indicating  a  thicker  boundary  layer  at  the  lower  tunnel  pressure. 

Temperature-sensitive  paint  (TSP)  runs  were  made  to  determine  the  undisturbed 
heating  rate  distributions  over  the  surface  of  the  booster.  Frames  from  motion  pictures 
taken  using  three  cameras  during  a  TSP  run  are  shown  in  Fig.  77.  For  this  run  (group 
194),  pQ  =  850  psia,  the  phase  change  paint  temperature  chosen  was  573°R.  This  relatively 
low  temperature  paint  melted  on  the  booster  nose,  canard  and  wing  leading  edges  before  the 
model  reached  the  tunnel  centerline.  Higher  heating  rates  are  also  evident  where  the  wing 
and  canard  shocks  strike  the  bod  /  and  in  the  wake  region  behind  the  canards.  The  span- 
wise  heating  rate  across  the  wi  jg  and  fuselage  is  seen  to  agree  qualitatively  with  the  span- 
wise  distribution  plotted  in  Fig.  74. 

Heat  transfer  and  pressure  distributions  were  obtained  along  rays  on  the  lower  sur¬ 
face  of  the  orbiter  (see  Fig.  16).  The  undisturbed  distributions  were  obtained  by  locating 
the  orbiter  high  above  the  booster,  above  the  booster  bow  shock,  and  measuring  the  heat 
transfer  rates  and  surface  pressures  for  orbiter  angles  of  attack,  (“  orjj )»  of  5  and  10 
degrees. 

Distributions  of  heating  rates  along  two  rays  on  the  lower  surface  of  the  orbiter,  for 
the  highest  and  lowest  tunnel  flow  stagnation  pressures,  are  indicated  in  Fig.  78  for 
a o^  =  °*  At  the  higher  tunnel  pressure,  the  inboard  ray  heating  rates  are  as  much  as 
four  times  as  large  as  the  heating  rates  measured  along  the  outboard  ray,  whereas,  at  the 
lower  tunnel  pressure,  the  heating  rates  along  the  inboard  ray  are  less  than  those  along  the 
outboard  ray.  The  pressure  distribution  along  the  centerline  of  the  lower  surface  of  the 
orbiter  is  indicated  in  the  last  portion  of  Fig.  78. 


I 


t; 


Corresponding  heat  transfer  and  pressure  distributions  for  Oq^  =  10°  are  indicated 
in  Fig.  79.  hi  these  cases,  the  outboard  ray  heating  rates  are  consistently  higher  than 
those  along  the  inboard  ray.  The  pressures  along  the  centerline  are  nearly  constant  for 


Curves  were  faired  through  the  measured  data  values  and  used  as  the  undisturbed 
reference  conditions  to  indicate  the  effects  of  interactions  caused  by  the  shock  generators 
on  the  booster  and  by  the  booster  on  the  orbiter. 

b.  Interaction  Effects  on  Booster 

Changes  in  the  heat  transfer  rate  distributions,  resulting  from  the  interaction  flow 
caused  by  the  0.5"  diameter  hemisphere-cylinder,  are  indicated  in  Fig.  80  for  a  tunnel 
flow  stagnation  pressure  of  850  psia.  Along  the  centerline  (row  1),  the  heating  rate  drops 
below  the  undisturbed  value  at  x  =  12  in. ,  and  then  rises  sharply  downstream  of  x  =  15  in. 
Along  row  2,  which  is  30°  off  the  centerplane  of  the  booster,  the  heating,  rate  starts  to 
drop  below  the  undisturbed  value  at  x  =  15  in. ,  and  doesn’t  rise  above  the  undisturbed 
value  until  x=sl7  in.  The  spanwise  distribution  of  heating  rates  is  affected  only  near  the 
centerplane  of  the  booster. 

The  orbiter,  which  has  a  0. 5  inch  diameter  spherical  nose,  affects  the  booster  heat¬ 
ing  distribution  much  more  extensively  than  the  smaller  hemispherical  shock  generator 
even  when  relatively  far  away  from  the  booster.  At  a  separation  distance,  Z^,  of  two 
inches,  the  booster  bow  shock  reflects  from  the  lower  surface  of  the  orbiter  and  interacts  1 
with  the  boundary  layer  on  the  booster  surface.  These  effects  are  evidenced  in  the  heating 
rate  distributions  shown  in  Fig.  81  and  also  in  the  corresponding  schlicren  flow  photographs. 
The  heating  rates  along  both  rows  1  and  2  drop  below  the  Undisturbed  values  at  k,  =  8  in.  The 
interference  heating  rates  are  considerably  larger  than  the  undisturbed  values  at  x  =  12  in. , 
and  exhibit  a  second  drop-then-rise  near  the  booster  base.  The  spanwise  distribution,  from  the 
fuselage  centerline  outboard  on  the  wing  surface ,  is  also  substantially  affected  by  the  interaction. 

<  i 

As  shown  in  Fig.  82,  the  interference  heating  is  greatly  aggravated  when  the  orbiter 
is  pitched  to  a  10°  angle  of  attack.  The  drop  in  heating  rates,  which  is  associated  with  the 
onset  of  separation,  occurs  upstream  of  x  =  8  in.  Both  streamwise  distributions  of  heating 
rates  attain  maximum  values  of  approximately  h™  =  0. 4  h„_  „  at  x  =  16  in.  These  values 

1  KIj  r  i 

represent  magnifications  of  approximately  18  times  the  undisturbed  heating  rate  for  row  2 

o  * 

(30  off  centerplane)  and  approximately  53  times  the  undisturbed  heating  rate  for  row  1  (in 

the  centerplane).  The  spanwise  heating  rate  distribution  at  station  x  =  15.6  inches  indicates 
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severe  heating  on  the  fuselage,  which  is  in  the  region  of  reattachment,  and  a  decrease  in 
heating  outboard  on  the  wing  surface,  which  is  in  a  separated  flow  region* 

Finally,  the  most  extensive  region  of  separated  flow  and  greatest  heating  was  caused 
by  the  10°  wedge,  as  indicated  in  Fig.  83.  The  boundary  layer  flow  separates  from  the 
booster  upper  surface  at  x  =  7  in.  There  are  two  peaks  in  the  downstream  heating  dis¬ 
tributions,  one  at  x  =  13, 5  in.  and  another  at  x  =  16  in.  The  spanwise  heating  distribution 
remains  high  over  the  entire  extent  of  the  instrumentation  on  the  upper  surfaces  of  the 
booster  and  wing.  Reynolds  number  effects  on  the  heating  distributions  are  also  indicated 
in  Fig.  83.  The  filled  symbols  represent  data  obtained  at  the  highest  tunnel  flow  stagna¬ 
tion  pressure  whereas  the  open  symbols  represent  the  corresponding  data  obtained  for  the 
lowest  tunnel  flow  pressure  level. 

Frames  from  motion  pictures  taken  during  the  temperature-sensitive  paint  runs  for 
the  10°  wedge  and  10°  orbiter  and  booster,  for  a  tunnel  stagnation  pressure  of  850  psia, 
are  shown  in  Figs.  84  and  85,  respectively.  A  paint  having  a  phase  change  temperature 
of  810°R  was  chosen  fur  the  wedge  generated  interaction.  The  times  of  the  sample  frames 
and  the  corresponding  heat  transfer  rates  along  the  boundaries  of  the  melted  paint  regions 
are  indicated  in  Fig.  84,  Sample  values,  obtained  from  these  frames,  are  plotted  in 
Fig.  83.  The  paint  and  thermocouple  results  are  seen  to  agree  quite  well  along  the  booster 
centerplane.  However,  the  paint  data  indicate  a  dip  in  the  spanwise  healing  rate  on  the 
wing  surface  that  is  not  indicated  by  the  thermocouple  data.  For  this  particular  run  (group 
196),  it  appears  that  too  heavy  a  coat  of  paint  was  sprayed  on  portions  of  the  booster 
surface,  which  may  have  caused  this  effect. 

Along  the  wing  leading  edge,  in  both  Figs.  34  and  85,  there  is  a  localized  inboard 
region  of  intense  heating  and. an  outboard  region  of  increased  heating.  The  inboard  "hot 

1 

spot"  is  attributed  hi.  self-induced  shock  impingement  resulting  from  boundary  layer 
separation  from  the  surface  ahead  of  the  wing-fuselage  junction  (Ref  17).  The  outboard 
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region  of  increased  heating  results  from  the  shock  generated  by  the  wedge  or  orbiter. 

A  paint  having  a  phase  change  temperature  of  960°R  was  chosen  for  the  10°  orbiter 
shock  interaction.  The  associated,  relatively  long  times  and  comparatively  high  heating 
rates  are  indicated  for  the  motion  picture  frames  shown  in  Fig.  85.  The  heating  rates 
from  these  frames  are  plotted  in  Fig.  82;  they  agree  well  with  tie  thermocouple  data. 

Changes  in  the  pressure  distributions  caused  by  the  10°  orbiter  and  10°  wedge  inter¬ 
actions  are  indicated  in  Fig.  86  for  the  highest  and  lowest  tunnel  flow  stagnation  pressure 

i 

levels.  At  the  higher  tunnel  pressure  level,  the  pressures  along  the  booster  centerline 

start  to  rise  to  a  plateau  level  at  x  =  7  inches  for  the  orbiter  generated  shock  interaction. 


The  pressure  peak  occurs  at  x  =  17  in.  Separation  occurs  somewhat  earlier  for  the  lower 
tunnel  pressure,  but  the  downstream  pressure  peak  still  occurs  at  x  =  17  in.  In  both 
cases,  for  the  orbiter  generated  shock  interactions,  the  centerline  pressure  distributions 
exhibit  "ell  defined  plateau  regions,  indicative  of  laminar  separation. 

Only  the  downstream  pressures  were  measured  for  the  wedge  generated  shock  inter¬ 
actions  shown  in  Fig.  86.  The  peak  pressures  are  larger  than  for  the  orbiter  generated 
shock  interactions,  and  occur  further  upstream.  As  uoted  in  the  figure,  the  pressure 
distributions  in  both  the  centerplane  (row  1)  and  in  a  plane  30°  off  centerplane  (row  3)  were 
essentially  identical  for  the  wedge  generated  shock  interactions. 

Streamwise  heating  rate  distributions  along  rows  1  and  2  on  the  booster  surface, 
measured  for  three  different  separation  distances  between  the  orbiter  (at a  =  0)  and  booster 
for  pQ  =  75  psia,  are  presented  in  Fig.  87.  As  the  separation  distance  decreases,  inter¬ 
action  effects  occur  closer  to  the  booster  nose,  and  there  are  multiple  peaks  in  the  heating 
distributions.  These  peaks  correspond  to  the  multiple  shock  reflections  and  interactions 
between  the  booster  and  the  orbiter.  Similar  peaks  are  evident  in  the  centerline  pressure 
distributions  shown  in  Fig.  88. 

The  half  inch  hemisphere-cylinder,  when  placed  near  the  booster  nose,  increased  the 
heating  along  the  booster  centerline  from  x  =  2  inches  to  x  =  18  inches,  as  shown  in  Fig.  89. 
The  increased  heating  is  associated  with  many  shock  reflections  between  the  booster  and 
hemisphere  cylinder.  However,  the  heating  was  not  increased  along  row  2  nor  on  the  wing 
surface  outboard  of  y  =  2  inches.  Downstream  of  x=«16  in. ,  the  data  can  be  affected  by 
shocks  from  the  model  support  structure,  and  therefore  may  not  be  valid.  Conversely,  the 
10°  wedge  increased  the  heating  and  pressure  both  on  the  booster  centerline  and  in  the  plane 
30°  off  centerplane,  as  shown  in  Fig.  90. 

c.  Interaction  Effects  on  Orbiter 

Heat  transfer  rates  and  pressures  were  measured  along  rays  on  the  lower  surface  of 
the  orbiter.  Changes  in  the  heating  rates,  caused  by  the  interaction  flow  between  the  orbiter 
(at  a=  0)  and  the  booster  at  three  separation  distances,  can  be  ascertained  by  comparing 
Figs.  78  and  91.  At  the  greatest  separation  distance,  the  heating  rates  dip  below  the 
undisturbed  values  on  the  orbiter  and  then  rise  to  peak  values  downstream  of  x  =  6  inches. 
The  amplification  is  most  severe  along  ray  4,  which  is  at  an  angle  of  4°  with  respect  to  the 
orbiter  centerline.  For  example,  at  x  =  7  inches  the  interference  heating  along  ray  4  is 
8  times  larger  than  the  undisturbed  value.  The  interference  heating  along  the  most  outboard 
ray  (No,  1),  is  50  to  75  percent  larger  than  the  undisturbed  heating.  The  locations  of  peak 
heating  move  forward  as  the  separation  distances  decrease.  At  the  closest  distance,  there 
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are  multiple  peaks  in  the  heating  distributions,  corresponding  to  multiple  shock  reflections 
between  the  booster  and  orbiter. 

Similar  to  the  undisturbed  flow  heating  rates,  the  interference  heating  rates  on  the 
orbiter  are  affected  strongly  by  changes  in  the  tunnel  flow  stagnation  pressure.  These 
effects  are  evident  by  comparison  of  Figs.  91a  and  92a  for  ^  =  0,  and  by  com¬ 
parison  of  parts  b  and  c  of  Fig.  92  for  «ork  =  10°.  The  heating  amplification  factors  are 
several  times  larger  for  the  higher  tunnel  flow  stagnation  pressures.  At  zero  angle  of 
attack,  the  amplification  of  the  heating  rate  along  the  centerline  is  doubled  and  the  amplifi¬ 
cation  of  the  heating  rate  along  ray  2  is  tripled  at  the  higher  tunnel  pressure  level.  At  a 
ten  degree  angle  of  attack,  the  peak  heating  rates  along  rays  2  and  3  are  4  to  5  times  larger 
for  pQ  =  850  psia  than  for  pQ  =  75  psia.  This  clearly  indicates  the  importance  of  Reynolds 
number  effects  for  shock  interactions  with  laminar  boundary  layers. 

As  expected,  the  «or^  =  10°  cases  result  in  higher  interference  heating  on  the  orbiter 
than  the  ^  =  0°  cases.  At  the  higher  tunnel  flow  pressure  level,  the  heating  along  both 
rays  2  and  3  was  essentially  identical. 

Data  obtained  from  the  corresponding  temperature-sensitive  paint  run  (Fig.  85), 
agree  well  with  the  thermocouple  data,  as  shewn  in  Fig.  92. 

Centerline  pressure  distributions,  corresponding  to  the  heat  transfer  rate  distribu¬ 
tions  shown  in  Figs.  91  and  92,  are  shown  in  Fig.  93.  The  multiple  shock  reflections  that 
occur  as  the  orbiter  and  booster  move  closer  together  are  evident  in  Fig.  93.  In  part  b, 
the  influence  of  tunnel  flow  pressure  level  is  indicated  clearly.  At  p  =75  psia  and  at 
pQ  =  850  psia,  the  pressure  peaks  at  approximately  x  =  6  inches.  In  both  cases  the  initial 
pressure  is  higher  than  the  undisturbed  value  (cf.  Fig.  79). 

4.  DATA  LOG 

Table  III  summarizes  the  significant  data  obtained  from  each  test  run.  The  table  is 
organized  by  type  of  receiver  (i.e. ,  sharp  flat  plate,  booster,  etc.)  and  then  type  of  gen¬ 
erator  (i.e. ,  wedge,  cone,  etc.).  Heat  transfer  data  obtained  from  thermocouple  and 
paint  runs,  and  pressure  data  are  presented. 

For  each  test  configuration,  the  table  identifies  the  te3t  group  numbers,  the  stagna¬ 
tion  pressure,  the  generator  and  the  nominal  location  for  shock  intersection  with  the  re¬ 
ceiver  (Nominal  x^).  In  addition,  the  location  at  which  the  impingement  shock  strikes  the 
boundary  layer  edge  (xBL),  and  the  actual  location  of  shock  intersection  with  the  receiver 
plate  (xp|atg)  are  presented.  These  dimensions  were  measured  from  the  profile  photographs. 
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Also  noted  for  each  run  is  the  location  of  peak  amplified  pressure  or  heating  rate,  (xpeaic)» 
and  the  undisturbed  pressure  or  heating  rate  at  xpea^.  Non-dimensional  pressures  were 
converted  to  dimensional  values  using  the  reference  properties  of  Table  IV. 


The  measured  heat  transfer  coefficients  based  on  gas  total  temperature,  hT,  were 

converted  to  the  more  commonly  used  coefficient  based  on  adiabadic  wall  temperature, 

h  ,  as  follows: 
aw 


h 


aw 


(T  -  T  ) 
v  o  w 

(T  -  T  > 
'  aw  w 


For  those  runs  where  the  undisturbed  boundary  layer  was  laminar  in  the  interaction 
region,  it  was  necessary  to  determine  whether  the  boundary  layer  was  laminar  or  turbulent 
after  the  interaction,  in  order  to  choose  the  correct  T  .  It  was  not  always  possible  to 
make  this  determination  by  examination  of  the  flow  photographs.  As  a  result,  plots  of 
amplification  in  hT  vs  amplification  in  pressure  were  constructed,  similar  to  those  dis¬ 
cussed  in  the  following  section,  and  these  were  used  to  discern  whether  the  disturbed 
boundary  layer  was  laminar  or  turbulent.  In  general,  where  the  data  plots  indicated  that 
the  disturbed  boundary  layer  was  transitional,  the  laminar  adiabatic  wall  temperature  was 
used  to  compute  h  .  In  those  cases  where  the  boundary  transition  was  judged  to  be  nearly 
complete,  the  turbulent  value  of  adiabatic  wall  temperature  was  used  to  compute  h  c  The 
subscript  L  or  T  with  the  value  of  h&w  denotes  whether  the  laminar  or  turbulent  adiabatic 
wall  temperature  was  used  to  compute  the  coefficient. 


For  some  runs,  the  pressure  or  heat  transfer  distribution  would  indicate  more  than 
one  apparent  peak  in  the  interaction  region.  In  these  cases,  the  first  peak  value  asso¬ 
ciated  with  the  primary  impingement  of  the  generator  shock  was  selected  by  comparison 
of  the  measured  distributions  with  profile  flow  photographs.  In  some  cases,  secondary 
peaks  could  be  correlated  with  reflected  shock  waves  or  with  the  generator  lip  shock  as 
evidenced  in  the  flow  photographs.  In  a  number  of  cases,  however,  the  selection  of  a  first 
peak  required  considerable  engineering  judgement.  In  these  cases,  the  secondary  peaks 
are  also  noted  in  Table  III. 


As  previously  mentioned  (Section  III— 1),  the  heat  transfer  data  were  obtained  at 
wall  temperatures  significantly  lower  than  for  the  pressure  data.  Usually,  this  disparity 
produced  no  clear  disagreement  in  flow  fields,  as  evidenced  by  the  profile  flow  photo¬ 
graphs,  pressure  and  heat  transfer  distributions.  In  some  cases,  however,  pressure  and 
heat  transfer  distributions  obtained  for  identical  test  configurations  showed  inconsistent 
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shapes,  number  or  location  of  peaks,  or  other  evidence  to  suggest  that  there  were  u^orent 
boundary  layer  flows  for  the  pressure  and  heat  transfer  runs.  Flow  photographs  (shadow¬ 
graphs  for  heat  transfer  runs  and  scblieren  photographs  for  pressure  runs)  for  these  runs 
were  also  compared.  These  runs  are  identified  by  appropriate  comments  in  Table  in. 

The  phase  change  paint  testing  provided  useful  qualitative  information  about  the  dis¬ 
tribution  of  heating  over  the  receiver  body.  This  type  of  testing  is,  however,  not  ideally 
suited  to  the  determination  of  peak  heating  values.  When  the  region  of  peak  heating  is 
small,  as  was  true  in  our  testing,  it  is  extremely  difficult  to  discern  the  exact  time  phase 
change  begins,  and  also  to  determine  the  location  of  the  point  of  peak  heating.  In  addition, 
paint  phase  change  temperatures  were  selected  to  provide  information  on  the  distribution 
over  the  receiver  body  in  a  reasonably  short  test  time.  Such  a  selection,  at  times,  re¬ 
sulted  in  the  occurrence  of  phase  change  at  the  location  of  peak  heating  shortly  after  the 
model  reached  the  tunnel  centerline.  This  introduces  an  error  in  evaluating  hT 
because  hT  is  not  constant  over  the  injection  period  although  it  is  assumed  constant  in  the 
solution  described  in  Section  II.  In  fact,  hT  varies  as  the  model  passes  through  the  tunnel 
boundary  layer.  As  this  transient  time  decreases  with  respect  to  the  time  required  for 
phase  change  to  occur,  the  error  decreases. 

As  a  consequence,  the  phase  change  paint  data  is  included  in  Table  III  to  indicate 
the  general  agreement  between  paint  data  and  thermocouple  data.  The  quantitative  data 
obtained  from  the  phase  change  paint  testing  have  not  been  used  in  the  analytical  correlations 
discussed  in  the  following  section. 
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SECTION  IV 
DATA  CORRELATION 

1.  ANALYTICAL  CONSIDERATIONS 

Although  it  is  not  our  intent  to  perform  rigorous  analyses  of  the  interaction  pro¬ 
cess,  some  discussion  of  a  very  simple  analytical  approach  provides  the  insight  nec¬ 
essary  to  organize  the  test  data  for  correlation.  The  analytical  assumption  we  will 
make  is  that  the  peak  heating  after  shock  impingement  can  be  represented  by  an  ideal 
zero  pressure  gradient  flat  plate  film  coefficient  equation.  This  may  be  valid  for 
either  of  two  reasons:  (1)  the  boundary  layer  may  actually  have  low  pressure  gradients 
in  the  x  and  y  directions  at  the  point  of  peak  heating,  since  this  is  some  distance  down¬ 
stream  of  the  impingement  point;  or  (2)  although  pressure  gradients  exist,  resolving 
the  boundary  layer  into  an  equivalent  flat  plate  flow  may  be  sufficiently  accurate  to  pre¬ 
dict  heating.  The  standard  forms  of  the  laminar  and  turbulent  flat  plate  heating  equations 


Laminar 


Turbulent 


Nu  =  .  332  Re* 5  Pr* 33 
Nu  =  .  0296  Re* 8  Pr* 33 


Simple  forms  of  these  equations,  from  handbooks  such  as  Ref.  18  ,  are 
u  _  .  00963  (pV)' 5 

"LAM"  ~  .tf.5  (8) 

-1*  X 

_  .  0334  (pV)*  8 

nTURB  t“.576  x  .2  (9) 

where  h  is  in  BTU/ft2  hr  °R,  p  is  in  LBS/FT2,  V  is  in  ft/sec,  T  is  in  °R,  and  x  is  in  feet. 
These  equations  result  from  expressing  density  by  the  perfect  gas  law  and  by  substituting 
functions  of  the  reference  temperature,  T*,  for  thermal  conductivity  and  viscosity.  The 
reference  temperature  is  defined  as 

[Tw+Te]+.22  [tr-T.] 

Using  subscript  1  for  the  initial  boundary  layer,  2  for  conditions  after  the  incident 
shock  and  3  for  conditions  after  the  reflected  shock,  we  can  write  three  possible  equations  for 
the  peak  heating  amplification  depending  on  type  of  boundary  layer. 

Case  1.  Boundary  layer  laminar  before  and  after  shock  impingement. 
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Casa  2.  Boundary  layer  turbulent  before  and  after  impingement. 


.0334  <p3  v3y 


3TURB_ 


m  .  576  .2 

x*  xo 

3  3 


*1  O 

TURB  .0334  (PjVj)  *8 


m  • 576  .  •  2 

x*  xi 


Case  3.  Boundary  layer  laminar  before  impingement  and  turbulent  after  impingement. 


.0334  (p3V3)- 
t  •  576 


*o 

°TURB_ 


*LAM  .00963  (pj  V2)  ,5 
Z  704  75 

T*i  *1 


=  3.47 


<Pi  V 

(T*  )  ’ 
1 


/p3y  /V 

l  pT)  IfT” 

V1/  \  a 


Our  testing  was  performed  at  M  =  8  with  a  range  of  shock  strengths  up  to  a  15°  wedge. 

Table  VII  shows  the  computation,  for  5,  10,  and  15  deg.  wedge  shocks  of  the  various  flow 

properties  required  in  equations  (10),  (11)  and  (12).  In  order  to  simplify  the  equations,  the 
T 

calculated  ratios  and  Vg  can  be  approximated,  by  a  curve  fit. 


as  powers  of  the  pressure  ratio: 


i/S 


vi  "i, 


r.  032 


-.  022 
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Substituting  these  ratios  in  equations (10),  (11,  and  (12)  and  using  our  Tunnel  B  values  for 

.5 


V-  and  T*  in  equation  (12)  we  get 

*  .48 

n„ 


I  AM 


LAM 


(13) 


TURB 

‘l 

1TUBB 

l3 

TURB 

li 

LAM 


.2 


.78 


1 


where  x  is  in  ft  and  p  in  lb3/ft 


(14) 


(15) 


In  all  these  equations,  x3,  the  distance  from  the  apparent  start  of  the  boundary  layer 
that  exists  after  impingement  must  be  determined.  Since  we  are  assuming  that  ideal  zero 
pressure  gradient  boundary  layers  exist  both  before  and  after  impingement,  jar  de¬ 
termine  the  boundary  layer  thickness  at  station  3  if  we  make  some  assumpi  so*>  regarding 
conservation  of  mass  and/or  momentum.  As  a  first  approach  the  simplest  sumption  is 
that  the  interaction  process  does  not  add  mass  to  the  boundary  layer,  and  that  the  mean 
boundary  layer  density  and  velocity  are  proportional  to  the  edge  velocity  and  density  for 


boundary  layers  of  the  same  type. 

piVi  = 

j*  v  4 

3V3  3 

(16) 

*1  P3V3 

P3V3T1 

(17) 

*3  P1V1 

P1V1T3 

The  razio  of  boundary  layer  thicknesses  before  and  after  impingement  can  be  evaluated 
from  the  flow  properties  in  Table  VII.  With  our  simplifying  assumption,  it  is  the  same  for 
either  laminar/laminar  or  turbulent/turbulent  interactions. 


*  1  =  2.87  5  deg  wedge 

*3 

=  5.87  10  deg  wedge  ’ 
=  8.09  15  deg  wedge 


for  M^ 


8 
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It  should  be  noted  that  v/e  do  not  claim  that  these  are  the  actual  ratios  of  boundary 
layer  thickness,  only  that  these  are  the  ratios  to  be  used  in  an  equivalent  ideal  flat  plate 
heat  transfer  equation  if  it  is  assumed  that  no  mass  is  added  during  the  interaction.  Mass 
addition  would  result  in  lower  values  of  4  /  6  . 

Xj 

The  ratios  of  ~  to  be  ujed  in  equations  (13)  &  (14)  can  be  determined  from  standard 
x3 

boundary  layer  equations  using  the  known  value  of  4^/  4^.  These  ratios  are  calculated  for 
the  pressure  ratios  corresponding  to  5,  10  and  15  deg  wedges  and  are  shown  in  Fig  94.  By 
curve  fitting  i  asa  function  of  p_/p.  (see  Figure  94)  equations  (13)  and  (14)  can  be  simpli- 

Xr.  o  l 

fied  to  J 

h  \  *73 

3LAM  ~/P3  j  (18) 

hl  \P1  / 

XLAM  \  1 / 


TURB 


TURB 


Since  our  assumption  of  constant  boundary  layer  mass  flow  results  in  the  maximum 
value  of  »  these  equations  should  be  somewhat  conservative.  The  available  turbulent 


data  support  this,  since  Sayano  (Ref  2)  and  Neumann  and  Burke  (Ref  19)  obtained  a  correla- 
h«  /p„\  *8 

tion  of  tr-  =  l — —I  .  Holden  (Ref  20)  and  Kains  and  Keyes  (Ref  21)  correlated  the  data 

hi  vv 


/p3V85 

('pT/ 


.  It  will  be  shown  later  that  eqn  (18)  is  also  slightly  conservative  compared 


to  our  laminar  data. 


The  case  of  an  initially  laminar  boundary  layer  that  is  turbulent  after  the  shock 
is  not  easy  to  analyze  even  with  our  very  simplified  analysis.  Natural  transition  has  been 
studied  for  years  and  is  still  not  well  understood,  and  shock  induced  transition  is  much 
more  complex.  In  our  test  program  we  varied  the  tunnel  pressure  level  and  x.,  the 
location  of  impingement  on  the  plate.  When  these  two  parameters  result  in  a  local  Reynolds 
number  less  than  natural  transition  the  question  arises  as  to  what  shock  strength  is  re¬ 
quired  to  cause  transition  to  fully  turbulent  flow.  The  more  stable  the  boundary  layer  the 
greater  is  the  shock  strength  (pressure  ratio)  required  to  cause  transition.  Since  boundary 
layer  stability  and  transition  are  primarily  a  function  of  Reynolds  number  and  Mach  number. 
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for  a  flat  plate  without  strong  heat  transfer  effects,  we  can  assume  that  the  transition 
pressure  ratio  is  an  inverse  function  of  Reynolds  number  at  our  constant  Mach  number,  or 


f  r" 


In  order  for  eqn  (20)  to  also  be  valid  for  natural  transition  it  must  satisfy  the  constraint  that 
as  Ret -♦  Re_,  p_/p  -♦  1.  Natural  Transition  data  on  a  plate  of  similar  bluntness, 


“'J>»  P^,/Pj 


tested atM  =  8  in  the  same  tunnel  is  presented  in  Ref  22  and  can  be  correlated  approximately 
6 

as  ReT  ~  p  .  In  Ref  23  this  unit  Reynolds  number  effect  was  attributed  to  aerodynamic 
noise  from  the  tunnel  wall  boundary  layer. 

Combining  the  natural  transition  constraint  and  eqn  (20)  we  get 

It  will  be  shown  later  that  a  functional  relation  of  this  form  provides  a  reasonable  correlation 
of  our  shock  induced  transition  data. 

We  have  now  analytically  developed  the  form  of  the  equations  to  predict  laminar  or 
turbulent  heating  amplifications  (eqns  18  &  19)  and  to  predict  the  pressure  ratio  required  to 
cause  transition  (eqn  21).  The  remaining  important  parameter  is  the  apparent  start  of 
the  turbulent  boundary  layer  after  transition  (xg)  which  is  required  in  eqn  (15)  for  laminar/ 
turbulent  interactions.  Unfortunately,  modifications  of  the  simple  flat  plate  approach  of 
equating  laminar  and  turbulent  momentum  thicknesses  at  the  transition  point  in  order  to 
solve  for  xg  did  not  agree  with  our  data  when  transition  was  shock  induced.  Therefore, 
later  in  this  section,  we  will  obtain  a  data  correlation  for  x„  but  an  analytical  rationale 
for  the  equation  requires  further  study. 

2.  SHARP  PLATE/WEDGE  TURBULENT  INTERACTIONS 

Figure  95  presents  the  peak  heating  as  a  function  of  the  measured  pressure  am¬ 
plification,  P„/p,  •  K  is  also  seen  that  the  measured  pressure  amplifications  agree  well 

O  * 

with  the  theoretical  amplifications  predicted  by  an  inviscid  calculation  across  the  incident 

*  8 

and  reflected  shock.  The  peak  hea'  be  predicted  by  h3  ^Pg\ 

\  ~  \?7// 

which  is  consistent  with  previous  expei-ments. 
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3.  SHARP  PLATE/WEDGE  INITIALLY  LAMINAR  INTERACTIONS 

Figure  96  shows  our  wedge  test  data  and  best  fit  correlations.  It  can  be  seen  that 

the  boundary  layer  stays  lairinar  for  shock  strengths  up  to  p  /p  ~  4  to  20  depending  on 

«  i  — 

initial  Refolds  number.  The  laminar  data  correlate  as/h^  \  /  pg \* 7  which  is  consistent 

VH/  W) 

with  our  simplified  analysis. 

The  transition  regime  between  laminar/laminar  and  laminar/turbulent  interactions  is 
defined  by  only  a  few  data  points,  so  our  correlations  are  indicative  of  the  trenc  but  certainly 
not  exact.  However,  using  our  correlation  from  Figure  96  to  define  the  shock  strength 
necessary  to  have  fully  turbulent  flow,  PT/Pj,  we  obtain  Figure  97.  Since  we  did  not 
obtain  natural  transition  on  our  plate,  two  data  points  from  Ref  22  are  shown  for  natural 
transition  at  pQ  =  850  psia  with  leading  edge  bluntness  of  .  007  in.  and  .  010  in.  Our  leading 
edge  is  .  007  in.  thick  but  is  flat-faced  instead  of  cylindrical,  thereby  probably  causing 
a  greater  bluntness  effect.  One  of  the  questions  raised  by  our  transition  data  is  that  the 
data  of  Gulbran  et  al  (Ref  24)  taken  in  the  same  tunnel,  was  found  by  Newman  and  Burke 
(Ref  19)  to  be  fully  turbulent  at  p./p.,  -  3  at  p  =  200  psia  and  xr»20  in. ,  compared  '  j 

u  1  O 

PT/pi  *8.5  predicted  from  our  data  in  Figure  97.  Their  plate  was  very  sharp,  however, 

with  natural  transition  occuring  at  x  os  3  ft  at  pQ  =  850  psia  based  on  the  data  of  Ref  22. 

Using  our  curve  fit  correlation  of  p^  ~  Pi  ‘ 8  we  can  see  from  Fig.  97  that  transi¬ 

re 

tion  at  x  a*  20  in.  on  a  very  sharp  plate  should  occur  for  PT  *  2. 8  at  p  =  200  psia. 

P1 

Therefore  we  can  say  that  our  transition  correlation  (eqn  22)  from  Fig.  97  is  reasonably 
consistent  with  shock-induced  data  from  two  experiments,  and  with  natural  transition  data  on 
plates  with  slight  leading  edge  bluntness  taken  in  the  same  wind  tunnel. 


1 


Prp 

P1 


„  .  2  .  5x, 
Pj  e  1 


(22) 


Also,  the  correlation  is  similar  to  the  form  that  we  expected  from  our  simplified 

5x  5 

analysis  since  over  the  .  5  to  2. 0  ft  range  of  x  in  our  test,  e  ’  ~  x  '  .  The  transition 

correlation  is  therefore  approximately 
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~  .2  .5 

P1  h 


or  PT 


which  agrees  with  eqn  (21). 

The  heating  rate  amplification  data,  for  a  boundary  layer  that  is  laminar  before  the 
shock  and  turbulent  after,  correlates  from  Fig.  96  as 


™  =  1.9  Pl 


This  correlation  can  be  compared  to  the  analytical  expression,  eqn  (15),  to  obtain  an  equation 
for  x3,  the  apparent  start  of  the  turbulent  boundary  layer,  however,  a  modification  of  eqn 
(15)  to  account  for  leading  edge  bluntness  should  be  made  first.  The  ideal,  zero  pressure 
gradient,  laminar  flat  plate  equation  used  to  derive  eqn  (15)  is  h^^  =  .46  Pl  *  3  *s 


obtained  from  eqn  (8)  by  substituting  our  tunnel  flow  conditions.  Figure  98  presents  our 

actual  laminar  baseline  data  with  no  shock  impingement,  which,  because  of  our  slight 

leading  edge  bluntness,  correlates  as 

.47 

h  a  5  Pi  (24) 

hLAM  *5-~r 
-59 


Modifying  eqn  (15)  by  substituting  eqn  (24)  for  h^^  we  obtain 

^TURB  =  1. 27  p  * 33  x.  *  59/Vy  76  1 

xlam  yi  j  3 

Equations  (23)  and  (25)  can  be  equated  and  solved  for  x,. 


=  .  67  p 


.11  .04 

1  X1 


An  alternate  method  of  correlating  x3,  which  eliminates  the  need  to  accurately  define 


is  to  solve  directly  for  x_  from  the  turbulent  heating  rate,  h 


3  TUEB. 


This  is  done 


in  Figure  99  where  the  turbulent  data  (initially  laminar)  is  plotted  in  the  form  .  634  p  ' 8  vs 

O 

h3 

•  •  2  p  8  *"* 

P]_»  since  X3  *  =  •  634  3  •  is  the  ideal  turbulent  flat  plate  equation,  eqn  9,  for  our  tunnel 

h3 

conditions.  The  correlation  obtained  in  this  manner  be  comes 
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=  .6  p1 


Since  this  i3  a  more  direct  solution  than  eqn  (26)  we  will  use  eqn  (27)  as  our  data  correlation 
for  the  apparent  start  of  the  turbulent  boundary  layer. 

4.  THREE-DIMENSIONAL  SHOCK  GENERATORS 

In  order  to  determine  the  effect  of  an  incident  three-dimensional  shock  on  our  basic 
wedge  correlation,  we  tested  a  limited  number  of  cone  and  sphere  snock  generators.  The 
intersection  of  these  shocks  with  a  plate  produces  a  flow  field  with  pressure  gradients  away 
from  the  peak  pressure  location  in  the  y  direction.  Flow  in  the  ±  y  direction  will  tend  to 
reduce  the  peak  pressure  below  that  obtained  with  a  two  dimensional  shock  of  the  same 
strength.  Although  the  pressure  will  be  lower,  the  heating  should  be  higher  (for  a  given 
pressure)  due  to  three-dimensional  thinning  of  the  boundary  layer.  Both  of  these  effects 
should  be  more  pronounced  for  sphere  generators  than  for  cones  since  the  spherical  shock 
strength  decreases  with  downstream  distance.  This  will  cause  higher  pressure  gradients  in 
the  y  direction  than  would  occur  for  an  equivalent  shock  strength  cone. 

a.  Turbulent  Data 

Since  the  initially  laminar  test  data  are  complicated  by  separation  shocks  and  boundary 
layer  transition,  it  is  desirable  to  look  at  the  turbulent  data  first  to  determine  3-D  effects. 
Figure  100  shows  the  cone  data  with  the  wedge  correlation  of  h  /p_\  *  for  comparison. 


^  _/!3\ 


It  can  be  seen  that  the  peak  pressure  is  less  than  the  theoretical  value  obtained  by  an  in- 
viscid  calculation  of  a  conical  shock  reflected  as  an  oblique  shock.  However,  as  expected, 
the  heating  correlation  is  about  25%  above  the  wedge  correlation  when  based  on  the  measured 
pressure  ratio.  An  interesting  point  about  the  turbulent  data  is  that  the  3-D  effects  compen- 

D  8 

sate,  in  that  the  heating  can  be  predicted  by  _‘_3  _  ■  *  3\  *  providing  the  theoretical 
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pressure  ratio  is  used. 
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As  expected,  the  sphere  data  (Fig.  101)  exhibits  stronger  3-D  effects  than  the  cone 
data.  The  pressure  amplification  is  considerably  below  the  theoretical  prediction  so  that 
using  _^3  _  /^3  \  ‘ 3 

hi  Wn 


Theory 
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would  result  in  predicted  heating  rates  high  by  as  much  as  a  factor  of  two.  On  the  other 
hand,  using  b'  me-asur*  d  pressure  ratio  will  under  predict  the  heating  by  a  factor  up  to 
40%.  In  order  'o  determine  a  correlating  equation,  it  is  reasonable  to  assume  that  spherical, 
-ihock  heating  amplification  will  approach  heating  due  to  a  wedge  shock  as  the  diameter  of  the 
sphere  increases  or  as  the  height  of  the  sphere  off  the  surface  decreases.  Figure  102 

g 

therefore  shows  the  value  of  (h0/h1 )/ (p„/p. )  *  versus  JZ  where  Z  is  the  height  of  the  sphere 
centerline  above  the  plate.  The  result  is  a  correlation  for  turbulent  flow  . 
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where  P3/P1  is  the  inviscid  calculation  of  the  incident  spherical  shock  reflected  as  an  oblique 

shock.  This  equation  should  of  course  be  used  with  caution  since  it  is  based  on  only  fjour 

data  points.  No,:e  that  the  equation  predicts  that  for  a  large  sphere  close  to  the  surface 

(Z  D/2)  the  maximum  heating  is  ^3  ^  .85  f^3\  *8 
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b.  Initially  Laminar  Data 

The  effect  of  conical  and  spherical  shocks  on  an  initially  laminar  boundary  layer  is 
considerably  more  complex  than  the  turbulent  case.  Figure  103  shows  all  the  sharp  plate 
cone  data  compared  to  the  correlation  previously  discussed  for  the  wedge  data.  The  heating 
amplification  data  are  plotted  against  the  measured  (not  the  theoretical)  pressure  amplifica¬ 
tion.  Although  the  data  are  limited  it  appears  that  transition  occurs  at  a  lower  pressure 
ratio  than  for  two  dimensional  shocks  and  that  the  turbulent  heating  level  h  ./h,  is  less 
sensitive  to  impingement  location.  Earlier  transition  might  be  expected  since  for  a  given 

pressure  ratio  the  incident  shock  strength  must  be  higher  for  a  3-D  shock.  The  apparent 

2  15  ■  1 

lesser  dependence  on  x  implies  that  the  correlation  xg  '  ~  Pj  '  Xj  determined  for  two- 
dimensional  shocks  is  different  for  3-D  shocks.  The  spherical  shock  data  shown  in  Fig.  104 
is  very  similar  to  the  conical  data.  Transition  occurs  ac  a  lower  pressure  ratio  and  the 
heating  seems  to  be  less  dependent  on  x  than  2-D  shock  daia. 
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5.  BLUNT  PLATE  RECEIVER 

The  blunt  plate  wedge  data  are  shown  in  Fig.  105  and  compared  to  the  correlation 
of  the  sharp  plate/wedge  data.  Although  the  data  are  not  extensive,  some  trends  are 
apparent.,  The  heating  amplification  at  22  inches  is  similar  but  somewhat  lower  than 
for  the  sharp  plate.  At  the  front  of  the  receiver  plate,  the  generator  locations  were 

•  t  « 

selected  to  provide  impingement  at  a  nominal  7  in.  The  lower  local  Mach 
number,  however,  resulted  in  greater  wedge  shock  angles  than  would  occur  at  M  =  8. 

In  addition, the  plate  bow  shock  interacted  with  the  wedge  on  some  of  the  runs,  as  can  be 
seen  in  Figures  67  and  71.  The  net  result  of  these  effects  is  that  peak  heating  occurred 
at  about  x  =  4  inches  for  the  forward  wedge  runs.  The  heating  amplification  at  x  =  4 
inches  is  higher  on  the  blunt  plate  than  at  x  =  5  inches  on  the  sharp  plate.  Also,  the 
pressure  ratio  required  to  cause  transition  is  lower  on  the  blunt  plate,  which  is  as  ex¬ 
pected,  since  natural  .transition  occurs  at  a  smaller  x  (see  Fig.  59).  The  blunt  plate 
heating  data  at  x  =  22  in.  is  slightly  lower  than  for  the  sharp  plate,  which  is  consistent 
with  pqn  (23)  for  the  following  rei.sons.  Although  p-j  on  the  blunt  plate  is  twice  the 
'sharp  plate'value,  the  average  impingement  location  for  the  blunt  cases  was  x  =  19  in. , 

22  55 

versus  x  =22  inches  for  the  sharp  plate,  so  that  the  term  pj*  Xj  *  is  approximately 

O 

the  same.  However,  the  constant  in  eqn.  (23)  is  proportional  to  Vj*  (see  eqn.  12) 

•  ip  .  536 

which  is  somewhat  smaller  for  the  blunt  plate  due  to  the  reduced  local  Mach  number. 

At  the  front  of  the  plate  the  smaller  constant  in  eqn.  (23)  is  more  than  offset  by  the  term 
Pp  *  xi  * V  since  the  local  blunt  plate  pressure  is  about  five  times  the  sharp  plate  value. 
Therefore,  the  blunt  plate  heating  amplification  data  at  x  =  4  in.  is  about  25%  greater  than 
the  sharp  plate  data  at  x  =  5  in.  While  the  blunt  plate  data  qualitatively  agree  with  the 
sharp  plate  equations  the  data  were  not  complete  enough  to  justify  a  detailed  correlation. 

6.  BOOSTER  RECEIVER 

The  undisturbed  flow  field  on  the  booster  has  some  of  the  characteristics  of  both 

the  sharp  and  blunt  plates.  As  shown  in  Fig.  76,  the  high  pressure,  low  Mach  number 

» 

flow  in  the  nose  region  expands  around  the  shoulder.  The  remainder  of  the  flow  over 
the  booster  is  at  nearly  the  free  stream  pressure,  with  a  low  pressure  gradient  which  is 
characteristic  of  a  sharp  flat  plate.  The  bow  shock,  however,  reduces  the  local  flow 
velocity  and  raises  the  static  temperature  which  results  in  lower  heating  at  a  given  x 
station  thau  occui,  on  the  sharp  plate.  The  booster  data  could  be  discussed  in  terms  of 
hg/h^  as  we  have  done  for  the  blunt  and  sharp  plates  by  deriving  an  expression  for  the 
unoisturbed  laminar  heating  rate  h.,  but  since  we  are  mainly  interested  in  the  shock 
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interaction  process  it  is  clearer  to  discuss  h,  directly.  The  first  step  is  to  determine  the 

O 

type  of  boundary  layer  present  downstream  of  the  shock.  The  10  deg.  wedge  was  tested  with 
the  booster  ac  pQ  =  850,  200  and  75  psia.  On  the  sharp  plate  (Fig.  96)  only  the  pQ  =  75  psia 
case  with  the  wedge  in  the  forward  location  was  laminar  after  a  10°  wedge  shock.  Since  the 
start  of  natural  transition  on  the  booster  at  pQ  =  850  psia  is  at  x  =  12  inches  (see  Fig.  74) 
and  on  the  sharp  plate  at  x>  24  inches,  we  can  assume  that  a  10°  wedge  shock  will  re¬ 
sult  in  a  turbulent  boundary  layer  for  all  booster  cases. 

As  discussed  in  Section  HI  the  theoretics’  inviscid  pressure  rise  for  a  10°  wedge 
shock  is  pQ/p-,  ~  20.  The  pressure  ratio  should  be  lower  on  the  booster  since  the  local 
Mach  number  is  less  than  8.  In  addition,  three  dimensional  effects  should  result  i"  a 
lower  pressure  ratio  for  impingement  on  the  forward  cylindrical  cross  section  (fuselage) 
than  on  the  aft  flat  plate  section  (wing)  of  the  booster.  The  data  (Table  III)  does  show 
these  effects  (neglecting  the  pQ  =  850  psia  nm  at  the  forward  location  which  has  an 
additional  pressure  amplification  due  to  the  reflected  booster  bow  shock). 

The  turbulent  heating  after  the  shock  can  be  predicted  from  the  standard  flat  plate 
equation  (9)  assuming  P3  is  known  from  test  data,  or  from  inviscid  shock  calculations. 
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h3  =  ‘°334  ZT576 
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(28) 
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The  ratio 


.576 


is  somewhat  lower  than  free  stream  conditions  due  both  to  the  lower 


local  Mach  number  on  the  booster  and  the  effect  of  the  10°  wedge  shock,  as  shown  in 
Table  VII.  Both  of  the  effects  are  relatively  small  so  we  will  use  free  stream  conditions 

.  .8 

(29) 


h3  =  .  634 


.2 


Although  natural  transition  occurs  at  a  smaller  x  on  the  booster  and  may  effect  the 
correlation  of  x^  we  will  apply  our  flat  plate  correlation  to  the  booster  runs. 


h3  =  1. 06 
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Using  the  data  from  Table  HI,  eqn  (30)  can  be  compared  to  the  measured  h^. 
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(30) 
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Wedge  Data/Booster  Receiver 


Nominal 

P1 

P3 

Data 

h3 

Eqn.  (30) 
h3 

x  Ft 

Lbs/Ft2 

Lbs/Ft2 

Btu/Hr-Ft2-°R 

Btu/Hr-Ft2-°R 

.5 

2.02 

23.6 

11.6 

12.9 

.5 

4.76 

58.4 

17.9 

23.4 

.  5 

17.1 

574.0 

106.0 

119.0 

1.16 

1.87 

26.6 

10.8 

13.1 

1.16 

4. 46 

63.4 

19.9 

22.9 

1.16 

16. 13 

200.0 

53.6 

47.7 

With  the  exception  of  the  last  data  point,  eqn.  (30)  overpredicts  the  data  by  10  -  30%. 


This  is  to  be  expected  since  the  value  of  _3_ 


.8 


used  in  eqn.  (30)  was  for  free  stream 

T  .576 
1  * 

conditions  and  is  lower  for  local  conditions  on  the  booster  due  to  the  bow  shock.  The 
last  data  point  is  at  pQ  =  850  psia  with  the  wedge  located  aft  on  the  booster,  so  that  the 
shock  impinges  at  x.  *  14  inches.  As  noted  above,  the  undisturbed  heating  is  transi¬ 
tional  at  this  location  on  the  booster.  This  may  explain  why  eqn.  (30)  underpredicts 
the  heating  since  the  correlation  of  xg  was  based  on  initially  laminar  data. 


7.  DESIGN  METHODS 


In  paragraph  IV-1,  we  postulated  relations  to  define  shock  induced  transition  and 
heating  amplification.  These  relations  were  not  based  on  rigorous  analysis  but  it  was 
hoped  they  would  indicate  the  parameters  of  importance  and  the  approximate  form  of  the 
equations.  In  paragraphs  IV-2  thru  IV-6  we  correlated  the  test  data  and  compared  the 
results  to  the  assumptions  of  paragraph  IV-1  in  order  to  establish  the  most  reasonable 
design  methods.  The  steps  recommended  to  obtain  an  engineering  estimate  to  an  inter¬ 
ference  heating  problem  are  as  follows: 

1)  Define  the  undisturbed  flow  field  including  local  velocity,  pressure  and 
temperature. 

2)  Specify  the  natural  transition  Reynolds  number 

3)  Determine  the  shock  strength  Pg/Pj  from  an  inviscid  calculation  or  test 
data  if  available 

4)  If  the  initial  boundary  layer  is  turbulent  h„  =  h. 

O  1 
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5)  If  the  initial  boundary  layer  is  laminar  estimate  the  pressure  ratio  required  to 
cause  transition  by 


.2  .  5x 

P1  e  1 


where  C  can  be  determined  from  the  natural 


transition  conditions,  =  pT>  =  xT  and  PT/p  =  1* 

pT  1 

6)  If  PQ/p,  <  .  5  —  assume  the  boundary  layer  after  the  shock  is  laminar  and 
«  1  P-j 

/p3\  -7 

predict  the  heating  by  h  =  h.  I  — ) 

J  XLAM  \pi/ 


PT 

7)  If  p„/p1  >  —  the  resultant  boundary  layer  is  assumed  to  be  fully  turbulent. 

O  1 

•  2  ^ ^ 

8)  Predict  the  apparent  start  of  the  turbulent  boundary  layer  by  x3  =  .  6  Pj  * 

V  *®  p  *® 

9)  The  turbulent  heating  is  then  given  by  h  =  0. 334  __3 _  ^3 

T*3  . 576  xg  . 2 

,  PT  pt 

10)  If  p  /p.  is  between  .  5  -r—  and  ~  the  boundary  layer  is  transitional  and  the 
3  1  P1  pi 

conservative  assumption  is  to  use  thp  turbulent  heating  rate  for  hg. 


The  above  approach  is  applicable  to  shocks  that  can  be  considered  two  dimen¬ 
sional.  Three  dimensional  shocks  appear  to  result  in  lower  heating  (if  based  on  the 
theoretical  2D  pressure  rise}  but  will  cause  an  initially  laminar  boundary  layer  to 
transition  to  turbulent  at  a  lower  shock  strength. 

It  should  be  noted  that  these  methods  are  based  on  correlating  data  at  M  =  8  and 

6 

a  Reynolds  number  range  of  .  4  to  3.7  x  10  /ft.  Caution  should  be  exercised  in  using  these 
equations  for  other  flow  conditions.  This  is  especially  true  of  the  expression  for  Xg,  for 
which  no  analytical  logic  has,  as  yeti  been  developed. 

8.  APPLICATION  OF  METHODS  TO  OR  BITER/BOOSTER 

We  obtained  data  on  several  very  interesting  orbiter/booster  configurations,  with 
complex  shock  patterns  impinging  and  reflecting  off  both  the  orbiter  and  booster,  that  can 
be  used  to  check  the  recommended  design  methods  in  the  previous  paragraphs.  Since 
this  data  was  obtained  very  late  in  the  program,  a  review  of  only  one  simple  case  can  be 
included  in  this  report,  although  all  the  data  is  presented  in  Table  in. 
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The  configuration  chosen  for  discussion  is  the  orbiter  located  aft  at  a  10  angle 
of  attack  to  the  booster.  The  tunnel  pressure  was  pQ  =  850  psia,  the  thermocouple  data 
are  from  group  403,  and  the  pressure  data  are  from  group  452.  The  data  are  presented 
in  Figures  82  and  86.  Since  we  are  primarily  concerned  with  the  analytical  prediction  of 
the  shoe1-  interaction  effects,  we  will  assume  that  the  undisturbed  flow  fields  o"  the 
orbiter  and  booster  are  known,  and  that  we  desire  a  conservative  preliminary  estimate 
of  the  heating  amplification  due  to  impingement  of  the  orbiter  shock  on  the  booster.  Since 
natural  transition  on  the  booster  begins  to  occur  at  x  «=»  12  inches  and  the  orbiter  shock 
for  this  configuration  will  impinge  at  x  =  16  inches,  we  can  assume  that  the  orbiter  shock 
is  strong  enough  to  separate  the  boundary  layer  and  cause  turbulent  heating  in  the  re¬ 
attachment  region.  The  undisturbed  booster  pressure  from  Table  III  is  pj  =  15. 8  lbs/ft2* 
As  discussed  in  Section  in,  Hill’s  correlation  of  the  laminar  separation  plateau  pressure 
is  Ppj  =  1.7  p^  =  26.9  Lbs/Ft2.  The  orbiter  shock  angle  is  13°,  for  the  10°  angle  of 
attack.  The  inviscid  pressure  rise  is  p^/P^j  «  11  at  M  =  8.  Although  the  local  Mach 
number  is  less  than  8  due  to  the  combined  effects  of  the  booster  bow  shock  and  the 
boundary  layer  separation  shock,  it  is  conservative  to  use  the  pressure  rise  for  M  =*  8. 
Therefore  pg  =  11  x  P  p  296  Lbs/Ft  .  Another  conservatism  in  the  prediction  of  the 
shock  pressure  rise  is  the  treatment  of  the  orbiter  shock  as  a  wedge  shock.  The  orbiter 
shock  will  be  somewhat  three  dimensional,  which  should  result  in  a  lower  pressure  rise. 
However,  in  the  case  of  cone  generators  we  have  shown  that  the  lower  pressure  rise  is 
offset  by  increased  heating  due  to  three  dimensional  thinning  of  the  boundary  layer  and 
the  theoretical  pressure  ratio  for  a  wedge  shock  should  be  used  to  predict  heating. 
Therefore,  using  p^  =  296  lbs/ft2,  p^  =  15.8  Lbs/Ft2  and  x  =  1.33  ft  in  eqn  (30) 
we  obtain 
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h3  =  1. 06 


.15  .1 


=  64 


Btu 


hr-ft2-°F 


which  compares  with  the  experimental  vau'e  of  h  =  60  Btu/hr-ft2-°F  from  Table  EH. 
Although  this  level  of  accuracy  is  probably  not  typical,  the  sample  calculation  does 
illustrate  the  steps  required  to  obtain  an  engineering  estimate  of  a  shock  interference 
heating  rate. 
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Section  V 

CONCLUSIONS 

Fairly  extensive  data  of  pressure  and  heating  amplification  due  to  wedge  shocks 
impinging  on  sharp  flat  plate  laminar  and  turbulent  boundary  layers  indicate  that  a 
simplified  analytical  approach  to  interaction  problems  is  feasible. 

1)  The  shock  induced  pressure  rise  can  be  predicted  by  an  inviscid  flow  field 
calculation  across  the  incident  and  reflected  shocks.  It  is  important  to 
include  the  possible  pressure  rise  due  to  flow  separation  and  to  consider 
extraneous  shocks  that  may  reflect  off  the  shock  generator  and  amplify  the 
pressure  rise  due  to  the  primary  shock.  A  fairly  common  cause  of  shock 
reflection  is  the  boundary  layer  separation  shock,  which  is  almost  always 
present  for  laminar  boundary  layers  and  for  turbulent  boundary  layer  inter¬ 
actions  with  strong  shocks. 

2)  Heat  transfer  for  laminar  interactions  can  be  predicted  from  the  pressure 


rise  by 


peak\_ 


/Ppeak  \ 
V  pUnd  / 


3)  Heat  transfer  for  turbulent  interactions  can  be  predicted  from  the  pressure 


rise  by 


peak\_ 


pUnd 


4)  The  prediction  of  heating  for  an  initially  laminar  boundary  layer  that  transi¬ 
tions  to  turbulent,  due  to  the  incident  shock,  requires  a  method  of  predicting 
the  apparent  start  of  the  resulting  turbulent  boundary  layer.  Correlations 
of  the  pressure  ratio  required  to  cause  transition  and  the  apparent  start  of 
the  turbulent  boundary  laTrer  are  presented  in  Section  IV  as  well  as  some 
logic  to  explain  the  form  of  the  pressure  ratio  correlation.  The  present 
simplified  analytical  approach  did  not  yield  an  explanation  for  the  correlation 
of  the  apparent  start  of  the  turbulent  boundary  layer.  This  theoretical 
understanding  is  necessary  before  the  correlation  can  be  applied  with  confidence 
to  conditions  other  than  those  tested. 
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Data  was  also  taken  with  conical  and  spherical  shock  generators  and  blunt  plate  and 
booster  receivers.  While  time  did  not  permit  a  detailed  correlation  of  this  data  some 
general  conclusions  were  reached. 

5)  Pressure  rises  due  to  3D  shocks  are  less  than  predicted  by  inviscid  calcu¬ 
lations  assuming  the  reflected  shock  is  two  dimensional.  The  heating  can 
be  conservatively  predicted  by  using  the  theoretical  pressure  ratio.  The 
3D  shock  appears  to  cause  an  initially  laminar  boundary  layer  to  transition 
at  a  lower  shock  strength  than  does  a  2D  shock. 

6)  The  heating  correlations  for  the  sharp  flat  plate  are  applicable  to  the  blunt 
plate  and  booster  providing  local  flow  conditions  are  used  where  values  are 
significantly  different  from  free  stream  conditions.  The  correlation  of  the 
pressure  ratio  required  to  cause  transition  must  be  adjusted  for  the 
different  natural  transition  points  of  the  blunt  plate  and  booster. 
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Table  I  Run  Schedule  Overview 


Entry 

Oate 

Occupancy 

Receiver 

Generators 

Type  of  Testing 

1 

11/30/71 

12/9/71 

43  hours 

Flat  Plate 
(Blunt  &  Sharp) 

Wedges,  Cones, 
Hemispheres, 

Orbiter  Vehicle 
(uninstrumented! 

1  hermocouple 

Booster  Vehicle 

Wedges,  Cones, 
Hemispheres, 

Orbiter  Vehicle 
(uninstrumented) 

Thermocouple 

Booster  Vehicle 

Wedges,  Cones, 
Hemispheres, 

Orbiter  Vehicle 
(uninstrumented) 

Phase  Change  Paint 

1A 

2/18/72 

8  hours 

Flat  Plate  (Sharp) 

Wedge,  Cones 

Thermocouple 

2 

2/24/72 

3/3/72 

56  hours 

Flat  Plate 
(Blunt  &  Sharp) 

Wedge,  Cones, 
Hemisphere, 

Orbiter  Vehicle 
(uninstrumented) 

Pressure 

Booster  Vehicle 

Wedge,  Cones, 
Hemisphere, 

Orbiter  Vehicle 
(uninstrumented) 

Pressure 

3 

3/29/72 

4/4/72 

32  hours 

Flat  Plate 
(Sharp) 

Wedge,  Cones, 
Hemispheres 

Phase  Change  Paint 

Booster  Vehicle 

Orbiter  Vehicle 
(instrumented) 

Thermocouple 

Booster  Vehicle 

Orbiter  Vehicle 
(instrumental) 

Pressure 

Booster  Vehicle 

Orb;ter  Vehicle 

Phase  Change  Paint 

Nominal  test  conditions  are  listed  below. 


Mach  Number:  8 


Tunnel  Total  Pressure,  Po  (PSIA) 

850 

400 

200 

75 

Tunnel  Tota1  Temperature,  T0(°R) 

1341 

1286  « 

1245 

1200 

Reynolds  number  per  foot: 

3.7x10® 

1.9x10® 

1.0x10® 

.4x10® 

T.  ble  til  Data  Summary  (Sheet  1  of  1 1) 
Sharp  ,-lat  PJ  ite  Initially  Laminar  BL  —  No  Trip 
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Table  III  Data  Summary  (Sheet  2  of  1 1) 
i.  Sharp  Flat  Plate  Initially  Laminar  BL  -  No  Trip 
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HEAT  TRANSFER  PRESSURE 


Table  III  Data  Summary  (Sheet  5  of  11) 
b.  Sharp  Flat  Plate  Initially  Turbulent  BL  -  Trip 
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Table  III  Oata  Summary  (Sheet  7  of  11) 


I  BTU  L  -  Laminar  ....  .... 

HR  FTZ'F  T  -  Turbulent  <49)  ®VB  •>  XPK’  h«w  PK 1  h»w.ljnit  *  163,01 

h**'Und  band  on  extrapolated  value  tor  a  lam  undid.  BL 
(47)  Tha  first  data  point  on  the  h-j*  plot  corresponds  to  at*PK 

(«)  C  up,  «nm  «« ttputs  P«k,  «h.  «tu.l  p«k  ,«V  to  mm,  ,W  ^pW^t.on  „  .nd.rtmpurfr.bl.  Horn  no» 

(46)  Group  360  is  a  repeat  of  359  with  the  gen  moved  aft  1.0"  (48)  Groups  345  &  346  art  repeats  Shock  is  too  weak  to  be  seen  on  shadowgraph 


Table  III  Oata  Summary  (Sheet  8  of  11) 
c.  Blunt  Flat  Plate  -  No  Tip 
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Table  III  Data  Summary 
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Table  IV  Refcrsnca  Conditions  for  Heat  Transfer  Data 


Stagnation  Pressure,  psia 

75 

200 

400 

850 

hpgp— Btu/hr— ft^*R 

31.5 

55.8 

78.4 

113.5 

M 

CO 

7.8 

7.92 

7.96 

8.01 

P  “  —  psfa 

1.21 

3.07 

5.93 

12.1 

To-»R 

1200. 

1245. 

1286. 

1341. 

TawL~°R 

1035. 

1073. 

1107. 

1154. 

TawT~OR 

1090. 

1130. 

1167. 

1216. 

T„-JR 

89.9 

92.0 

94.0 

96.9 

T  —  °R 
w 

535. 

535. 

535. 

535. 

T.  —  °R 

LAM 

520. 

529. 

537. 

549. 

T.  -°R 

TURB 

533. 

542. 

551. 

562. 

(To  -V/,TawL-Tw» 

1.33 

1.32 

1.31 

1.30 

ITo-Twl/ITawrV 

1.2 

1.19 

1.19 

1.18 

Re/ft 

4.09  +  05 

1.01  +  06 

1.88  +  06 

3.70  +  06 

-  ft/sec 

3653. 

3722. 

3772. 

3864. 

p„  -  lb  /ft3 

00  m 

8.10  -  06 

2.00  -  05 

3.77  -  05 

7.47  -  05 

p»  ”  lb/ft — sec 

7.24  -  08 

7.40  -  08 

7.57  -  08 

7.80  -  08 

Table  V 

Data  Obtained  on  Flat  Plate  with  Blunt  Leading  Edge 


Wedge  Angles 

Cone  Haif  Angles 

Hemispheres 

Orbiter  a 

1 

po 

X. 

1 

5° 

10° 

15° 

2.5° 

7.5° 

0.5" 

2.0" 

0° 

5° 

P  T 

P  T 

P  T 

T 

P  T 

P  T 

_ Lj 

P  T 

T 

850 

7 

• 

•  • 

• 

• 

•  • 

n 

•  # 

22 

•  • 

J 

L*  • 

M 

■D 

H 

1, 

400 

7. 

•  • 

■  ■ 

■1 

■ 

• 

• 

22 

•  • 

M 

1 

19 

• 

75 

7 

■ 

•  • 

3 

ns 

m 

• 

• 

mil 

22 

■H 

ns 

■ 

HI 

1 

■  I 

*  Data  obtained  for  two  or  more  heights  of  shock  generator  above  plate  surface. 


p  —  Pressure  T  -  thermocouple 


Table  VI 

Data  Obtained  on  Booster  and  Booster-Orbiter  Configurations 


Po 


GENERATORS 

10°  Wedge  |  7.5°  Cone  |  0.5"  Hemisph 


p  T  P 


P  T  P 


T  P 


a  =  0 
P  T  P 


ORBITER 
a  =  5° 

P  T  P 


a  =  10° 
p  T  P 


850 


FORE 

AFT 


•  •  • 
•  •  • 


•  • 
•  • 


•  • 

•  •  • 


200 


FORE 

AFT 


75 


FORE 

AFT 


•  • 
•  « 


*  Data  obtained  for  two  height,',  of  orbiter  above  booster  surface, 
p  —  Pressure  T  -  Thermocouple  P  —  Paint 
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Fig.  11  Orbiter  Generator  Installed  with  Booster  Receiver 
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b.  TYPICAL  SECTION -PRESSURE  PANEL 


Fig.  12  Typical  Thin  Skin  and  Pressure  Insert  Cross-Sections 
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Fig.  16  Orbiter  Mod-il  Pressure  Tap  and  Thermocouple  Locations 


PRESSURE  AND  THERMOCOUPLE 


Fi<i.  17  Booster  Model  Pressure  Tap  and  Thermocouple  Locations 


-~r-  -  i^v  -v^' '  T--^ 

— --J >**  'V  T  ,.V'''’*-?.’--:;  ‘  . v.-v-; :  .<*  *  ■»,'.  . - ,  !  ... 


3 i5By.”5?rTv^'xn’ci^?^i.7si 


CAMERA  0 
(LEFTSIDE  VIEW) 


a.  CAMERA  SET-UP  FOR  TEFLON  RECEIVERS  (FLA  PLATE,  BOOSTER)  AND  STEEL 
GENERATORS.  (GROUPS  194-202,  375-387) 


b.  CAMERA  SET-UP  FOR  TEFLON  BOOSTER  WITH  TEFLON  ORBITER  (GROUPS  454-472) 


Fig.  18  Tunnel  Cross-Section  Showing  Camera  Set-Up  for  Phaae  Change  Paint  Tests 
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b.  SCHLIERENS 


Fig.  26  Shadowgraph  and  Schliaren  Photographs  of  Undisturbed  Flows  on  Sharp  LE  Plate; 
p0  =  850  psia  (Groups  108  &  254) 
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b.  SCHLIERF.NS 


Fig.  32  Shadowgraph  and  Schlieren  Photographs  of  Interaction  Caused  by  15°  Wedge  Shock; 
pQ  =  75  psia,  Xj  =  5  in.  (Groups  96  &  317) 
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UNDISTURBED  PRESSURE 


a.  SHADOWGRAPHS 


b.  SCHLIERENS 


Fig.  36  Shadowgraphs  and  Schlieren  K  ' 
p0  =  850  psia,  Xj  =  22  in.  (Groups 


•  *  Interaction  Caused  by  10°  Wedge  Shock; 
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Fig.  38  Shadowgraph  and  Schlieren  Photographs  of  Interaction  Caused  By  15°  Wedge  Shock; 
pQ  =  850  psia,  Xj  =  22  in.  (Groups  218  8<  275) 
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Fig.  39  Boundary  Layer  Correction  for  Wedge  and  Cone  Generator  Shock 
Strength  Calculation* 
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CALCULATED  PRESSURE  AMPLIFICATION  (P-j/P^ 


Fig.  43  Pressure  Amplif;  ;tion  for  Cone  and  Hemisphere-Cylinder  Generators  Sharp  Flat  Plate  (No  Trip) 
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CALCULATED  PR  :SSURE  AMPLIFICATION.  P3/P, 


Comparison  of  Measured  Pressure  Ratios  for  Wedge  Generators  with  Tripped  and  Untripped  Boundary 
Layers  at  p„  «  850  psia  and  x:  <*  22  in. 
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b.  PRESSURE 


Fig.  53  Heating  and  Pressure  Distributions,  on  Plate  with  BL  Trip,  Caused  by  Raised  15' 
Wedge  Shock;  pQ  =  850  psia,  Xj  ■  22  in.  (Groups  210  &  259) 
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Comparisons  of  Heating  Distributions  on  Sharp  L.E.  Plate 
Caused  by  0.5  in.  Hemisphere  and  7.5°  Cone  for 
Turbulent  and  Laminar  Boundary  Layers 


a.  SHADOWGRAPHS 


h.  SCHLIEREN 


Fig.  GO  Shadowgraph  and  Schlioren  Photographs  of  Undisturbed  Flows  on  Blunt  LE  Plate; 
P0“  850  psia  (Groups  104,  105,  and  348) 
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a.  HEAT  TRANSFER 
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Fig.  66  Heating  and  Pressure  Distributions  on  Blunt  Plate  caused  by  10°  Wedge  Shock; 
pQ  *  850  psia,  Xj  *  7  in.  (Groups  46, 345  8t  346) 
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a.  HEAT  TRANSFER 


x,  IN. 

b.  PRESSURE 


Fig.  70  Halting  and  Pressure  Distribution*  on  Blunt  Plata  causad  by  1S°  Wadga  Shock;  p„  -  75 
psia,  X| a  7  in.  (Groups  136  &  359)  ° 
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Fig.  76  Undisturbed  Pressure  Distributions  on  Booste*’, 
P0  =  75  snd  850  psia  (Groups  363  &  372) 
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Fig.  80  Interference  Heating  Distributions  on  Booster  Caused  by  0.5  inch 
Hemisphere,  p  =  850  psia  (Group  157} 


Fig.  81  Interference  Heating  Distributions  on  Booster  Caused  by  Orbiter 
at  Oqrq  =  0,  p0  =  850  psia  (Groups  391  &  392) 
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Fig.  85  Frame*  from  Motion  Picture*  Taken  During  Temperature  -  Sensitive  Paim  Run.  Phaw 
Change  Temperature  =  960°R,  Booster  and  10°  Orbiter,  pQ  =  850  p*ia  (Group  462) 
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Fig.  88  Interference  Pressure  Distributions  on  Booster  Caused  by  Orbiter 
(°OR3  “  at  Two  Different  Separation  Distances  from  Booster, 
pQ  3  75  psia  (Groups  44?  and  444) 
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Fig.  91  Interference  Heating  Distributions  on  Orbiter  (°qp  -j  =  0)  Caused 
by  Booster  at  Three  Different  Separation  Distances  Between 
Orbiter  and  Booster,  p  =  75  psia  (Groups  411, 414, 417,  and  434-436) 
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Fig.  93  Interference  Pressure  Distributions  on  Orbiter  Centerline  Caused  by 

Boo«ter  (Oq^q  =  0  and  10°),  pQ  =  75  and  850  psia  (Groups  442-444, 

447,  and  452) 


Fig.  95  Sharp  Plate-Wedge  Data  Turbulent  Boundary  Layer 
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